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ABSTRACT 


Microseisms  at  five  observatories  viz.  Victoria, 
Penticton,  Edmonton,  Resolute  and  Ottawa  are  studied  with 
special  reference  to  their  content  of  physical  wave  types, 
mechanism  of  their  propagation  in  the  crust  and  the  direc¬ 
tion  of  approach. 

Correction  for  the  instrument  response  is  made  us¬ 
ing  a  transfer  function  of  fourth  order  in  frequency.  Plot 
of  this  function  is  matched  with  the  calibration  curves  sup¬ 
plied  by  the  observatories.  The  matching  shows  that  in  the 
case  of  Ottawa  only  the  critical  response  conditions  are  met; 
but  at  Victoria,  Penticton  and  Resolute  the  seismometer  is 
overdamped  and  the  galvanometer  underdamped,  while  at  Edmon¬ 
ton  the  seismometer  is  underdamped  and  the  galvanometer  over¬ 
damped  . 

Peaks  in  the  autopower  and  the  crosspower  spectra 
occur  at  around  6.5  seconds,  4.5  seconds  and  2.5  seconds. 

The  main  peak  is  at  6.5  seconds  at  all  stations  although  the 
sedimentary  layers  at  Edmonton  seem  to  have  a  pronounced 
effect  in  accentuating  the  shorter  period  peak  at  2.5  seconds. 
Coherency  between  the  vertical  and  horizontal  components  is 
in  general  very  low  (often  much  less  than  0.7)  while  that  be¬ 
tween  the  horizontal  components  is  usually  around  0.8  or 
greater.  The  phase  values  between  the  vertical  and  horizon¬ 
tal  components  strongly  suggest  that  body  waves  (perhaps  LP 
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modes  of  Phinney)  make  up  a  considerable  portion  of  micro¬ 
seisms.  This  is  reinforced  by  the  sequential  peaks  in  the  ve¬ 
locity  sensitivity  spectra  and  by  the  studies  on  synthetic 
microseisms.  White  noise  with  a  peak  to  peak  amplitude  twen¬ 
ty  times  that  of  the  signal  is  required  to  change  the  phase 
and  coherency  for  a  simple  Rayleigh  type  of  motion  by  the 
amount  observed  in  the  data.  The  same  amount  of  white  noise 
would  also  distort  the  shape  of  the  microseisms  as  recorded 
on  the  instruments.  But  since  the  microseisms  appear  as  sinu¬ 
soids  on  the  records,  it  is  felt  that  "noise"  associated  with 
Rayleigh  or  Love  type  of  motion  is  not  background  white  noise 
but  some  seismic  wave  with  the  properties  of  body  waves  and 
pseudo-surface  waves.  Interference  of  Rayleigh  waves  and  the 
long  period  body  waves  would  then  be  explained  as  a  reason 
for  low  coherency  between  the  vertical  and  horizontal  compon¬ 
ents  . 

Studies  on  directional  spectra  indicate  that  the  un¬ 
certainty  of  direction  of  approach  increases  as  the  distance 
from  the  coast  increases.  This  is  attributed  to  the  coastal 
effect  rather  than  due  to  the  source  of  microseisms.  It  is 
suggested  that  the  elastic  vibrations  set  up  in  the  ocean 
bottom  by  the  second  order  fluctuations  are  channeled  in  the 
oceanic  crust  and  travel  in  different  directions  until  they 
reach  the  coast  from  whence  the  waves  travel  in  the  continent 
along  different  paths  so  that  the  sites  near  the  coast  have  a 
smaller  angle  of  vision  of  the  source  of  elastic  energy  and 
the  farther  sites  have  a  larger  angle  of  vision  resulting  in 

greater  uncertainty. 
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CHAPTER  I 


INTRODUCTION 

1.1  Bertelli  (Macelwane,  1952)  seems  to  be  the  first  to 
have  systematically  studied  microseismic  phenomenon.  He 
did  this  by  means  of  several  thousands  of  observations  on 
free  pendulums  and  on  reflections  from  the  surface  of 
mercury.  Darwin  &  Darwin  (Macelwane,  1952)  had  to  abandon 
their  study  of  lunar  tides  in  the  solid  earth  because  of 
this  phenomenon,  which  interfered  with  their  experiments. 

Milne  (Macelwane,  1952),  on  the  basis  of  work  extended  over 
several  years  (1880-1887)  and  done  in  different  parts  of  the 
earth  (Japan,  Italy,  France  and  England),  concluded  that  the 
whole  globe  is  affected  by  these  phenomena  and  he  called 
them  "microseismical  disturbances."  He  was  also  the  first  to 
recognize  a  "microseismic  storm"  and  to  coin  the  phrase,  in 
addition  to  observing  that  these  disturbances  can  travel  for 
long  distances.  Meanwhile  E.  Von  Rebeur-Pachwitz  (Macelwane, 
1952)  on  the  basis  of  work  done  in  Germany  during  the  period 
1892-1894,  submitted  a  report  on  Earth  tides  in  1895  wherein 
a  section  was  devoted  to  observations  on  microseisms.  Omori 
(1901)  recognized  the  periods  of  these  pulsations  to  be  about 
4  sec.  and  8  sec.  By  this  time  interest  in  microseisms  was 
widespread  and  at  meetings  of  the  International  Seismological 
Association  (Macelwane , 1952 )  a  standing  committee  was  appointed 
for  the  study  of  microseisms. 
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1.2  Many  theories  have  been  proposed  to  explain  the  origin 
of  microseisms  but  the  surf  theory  seems  to  be  the  one  with 
some  historical  importance.  Weichert  (Macelwane,  1952)  pro¬ 
posed  this  theory  during  the  second  international  conference 
in  Strasbourg.  Surf  approaching  a  continental  boundary  im¬ 
parts  energy  to  it  which  is  then  transmitted  as  microseisms 
in  the  continent.  Gutenberg  (Macelwane,  1952),  in  a  doctoral 
dissertation,  vigorously  supported  the  surf  theory  on  the 

basis  of  a  strong  correlation  between  the  occurrence  of  micro- 

/ 

seisms  and  waves  approaching  the  Norwegian  Coast.  However, 
Banerji  (1930)  and  Gherzi  (Macelwane,  1952)  found  that  the 
strongest  microseisms  occurred  some  hours  before  the  storm 
reached  the  Coast.  To  meet  this  objection  the  surf  theory 
was  modified  to  state  that  the  wave  energy  produced  by  the 
storm  radiates  in  the  ocean  to  form  a  swell  until  it  reaches 
the  continental  boundary  where,  by  a  mechanism  not  speci¬ 
fied  in  the  theory,  it  is  imparted  to  the  continents. 

1.3  Imbo  (Macelwane,  1952),  by  comparing  waves  in  the  Medit- 
erranian  with  microseisms  in  Sicily,  found  that  the  period  of 
the  microseisms  was  half  that  of  the  ocean  waves.  A  similar 
result  was  reported  by  Deacon  (1947)  and  Darbyshire  (1950). 

This  clearly  can  not  be  explained  by  the  surf  theory. 

The  theory  of  oceanic  storms  as  proposed  by  Banerji 
(1929,  1930,  1935)  and  Gherzi  (Macelwane,  1952)  suffered 
from  the  objection  that  the  internal  pressures 


'  , 


3 

due  to  water  waves  in  deep  water  die  off  too  rapidly  with 
depth  to  be  effective  in  imparting  energy  to  the  ocean  floor. 

As  early  as  1831  Faraday  (Longuet-Higgins ,  1952) 
knew  that  fluid  resting  on  a  vibrating  elastic  plate  will 
form  itself  into  standing  waves.  He  also  showed  that  the 
period  of  the  standing  waves  in  the  fluid  was  twice  that  of 
the  elastic  vibrations  in  the  plate.  Rayleigh  (Longuet- 
Higgins,  1952)  repeated  Faraday's  experiment  in  1883  and  ver¬ 
ified  by  a  slightly  different  method  the  doubling  of  the 
period.  Neither  Faraday  nor  Rayleigh  evaluated  the  second 
order  terms  involved  in  the  analysis  of  standing  waves.  Miche 
(Longuet-Higgins,  1952)  using  a  Lagrangian  System  of  coordin¬ 
ates  noticed  the  unattenuated  terms  and  evaluated  them. 
Longuet-Higgins  (1950)  expanded  the  standing  wave  theory 
where  the  second  order  terms  do  not  vanish  with  depth  and  if 
the  water  area  concerned  is  large  they  could  cause  microseisms 
on  the  continents.  Theoretically  then  the  period  of  ocean 
waves  would  be  twice  that  of  the  microseisms. 

1.4  Longuet-Higgins  (ibid)  suggests  that  interference  of 

waves  travelling  in  opposite  directions  may  be  necessary  to 
cause  the  second  order  fluctuations.  The  opposite  travel¬ 
ling  groups  of  waves  may  occur  in  cyclones  or  by  reflection 
of  sea  waves  from  a  coast.  The  interference  of  these  waves 
should  set  up  a  standing  wave,  where  as  shown  by  Miche,  the 
second  order  fluctuations  do  not  attenuate  with  depth 
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and  are  independent  of  the  depth  of  the  liquid  body  con¬ 
cerned.  If  the  period  of  the  seismic  wave  is  T  and  its 
wavelength  L,  then  the  pressure  applied  on  an  area  suf¬ 
ficiently  small  compared  with  L  may  be  considered  to  have 
been  applied  at  a  point.  The  time  difference  involved  in 
applying  any  pressure  at  another  point  of  the  area  would  be 
small  compared  with  T.  The  seismic  wavelength  is  very  large 
compared  with  the  length  of  a  gravity  wave  (sea  wave)  of 
comparable  period.  Hence  it  is  appropriate  to  consider  the 
properties  of  the  mean  pressure  over  a  large  number  of 
wavelengths  in  any  gravity  wave. 

It  can  be  easily  shown  that  the  pressure  fluctuations 
in  a  progressive  wave  over  a  multiple  of  wavelengths  is 
too  insignificant  to  cause  seismic  waves  in  the  sea  bed. 

In  a  progressive  wave  the  contributions  to  the  disturbance 
from  different  parts  of  the  sea  bed  tend  to  cancel  one 
another  out.  A  second  reason  why  pressure  fluctuations  in 
a  progressive  wave  are  ineffectual  in  causing  seismic 
oscillations  on  the  sea  bed  is  that  not  only  the  mean  pres¬ 
sure  but  also  the  pressure  difference  decreases  rapidly 
(exponentially)  below  a  wavelength. 

Now  consider  two  progressive  waves  of  equal  wave¬ 
length  and  amplitude  travelling  in  opposite  directions. 

The  mean  pressure  on  the  bottom  between  two  nodal  points 
vanishes  to  a  first  approximation.  But  the  addition  is  not 
exact  and  the  superposability  is  not  exact  as  the  equations 
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of  motion  are  nonlinear.  Consider  the  mass  of  water  contained 
between  the  standing  wave  on  the  free  surface,  the  sea  bottom, 
the  two  anti-nodal  planes  and  a  unit  length  in  a  direction 
orthogonal  to  the  direction  of  wave  propagation.  Since  there 
is  no  flow  of  mass  across  the  nodal  planes,  this  mass  consists 
always  of  the  same  particles;  therefore  the  motion  of  center 
of  gravity  of  this  mass  is  that  due  to  external  forces  alone 
which  act  upon  it.  Comparing  the  motion  of  water  with  that 
of  a  pendulum,  it  can  be  seen  that  when  the  pendulum  passes 
through  the  center  of  oscillation  the  equivalent  surface  of 
the  water  is  relatively  flat  and  when  the  pendulum  reaches 
the  point  of  its  maximum  displacement,  there  will  be  a  cor¬ 
responding  crest  on  the  water  surface.  The  center  of  grav¬ 
ity  of  the  mass  is  higher  when  crests  are  formed  than  when 
it  is  relatively  flat.  The  center  of  gravity  is  thus  raised 
and  lowered  in  one  complete  cycle.  The  external  forces  act¬ 
ing  on  this  mass  of  water  are  1)  that  due  to  gravity  which 
is  constant,  2)  the  atmospheric  pressure  P  on  the  surface 
is  constant  or  it  gives  a  constant  downward  force  of  XP , 
where  X  is  the  wavelength,  3)  the  forces  across  the  vertical 
planes  which  must  have  zero  vertical  components,  the  motion 
being  symmetrical  about  these  planes,  4)  the  force  on  the 
bottom  which  equals  XP.  Since  all  the  forces  other  than  XP 
are  constant,  XP  must  fluctuate  with  time.  In  other  words 
p  must  fluctuate  with  time.  The  mean  pressure  fluctuation  is 
proportional  to  the  square  of  the  wave  amplitude  on  the  free 
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If  F  is  the  force  acting  vertically,  z  is  the 
vertical  coordinate  dependent  on  time,  m  is  the  mass  of 
the  water, 

2 

„  ym3  Z'  3  2  v 

F  =  E -  =  - Emz 

3t2  8 12 

In  an  incompressible  fluid 

Emz  =  pj"^2C2dx  +  constant 
o 

where  £  is  the  vertical  displacement  of  the  free  surface, 
x  is  the  horizontal  coordinate,  p  density  and  ^  the  wave¬ 
length.  But  the  force  F  =  A  (P-Po-  pgh) ,  p  is  the  den¬ 
sity,  P  the  mean  pressure  at  the  bottom,  h  is  the  mean 
depth  of  water,  Po  pressure  at  the  free  surface. 

The  standing  wave  can  be  represented  by 
£  =  a  Cos  kx  Cosat 

where  a  is  the  amplitude  of  the  wave,  k  =  2tt /\,o  =2tt/t, 

(t ,  being  the  period  of  the  wave) ,  higher  order  terms  being 
omitted. 

After  simplification  we  find 
(P~po)/p-gh  =  -1/2  a2a2Cos2at 
This  shows  that,  to  the  second  order,  the  mean  pressure  P 
fluctuates  sinusoidally  with  twice  the  frequency  of  the 
original  wave,  and  with  an  amplitude  proportional  to  the 
square  of  the  wave  amplitude.  The  pressure  fluctuation  is 
independent  of  depth,  for  a  given  frequency  of  the  wave, 
though  the  frequency  itself  is  dependent  on  the  depth, 


- 
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given  by 


o2  =  gk  tanh(kh) 

For  two  progressive  waves  we  have 

p-P0  -  gh  =  -2a  a  aCos2at 
~~z —  .  1  2 


For  an  incompressible  liquid 


P~P o  -  gh  =  -1/2  a2g.2Cos  (2a(z-h)/cl)  Cos  2at 

Cos (2ah/c ' ) 

r' 

very  nearly,  where  c'  is  the  velocity  of  sound  in  water. 
For  sea  waves. 


C(*/Y #t) 


,00  .oo 

Re  f  /  A(u,v)e' 

—  oo_oo  \  r  / 


(ukx+vky+ot)du  dv 


where  x,  y  are  the  horizontal  coordinates,  k  is  a  constant, 
and  a  is  a(u, v) . 


a2  =  gk/(u2+  v2)  tanh.  kh/(u2  +  v2) 

A  (u,v)  is  in  general  complex  and  Re  denotes  the  real  part. 

Longuet-Higgins  goes  on  to  prove  that  mean  pressure  vari¬ 
ation  for  a  three  dimensional  case  is 


P-Pq 

P 


l/2gX 

g 


9  2 
3t* 


1 

4R7 


-J.J 


00 


1/2  C 


2dxdy 


where 


=  Re  rr  A'  (u,v)ei(ukx+vky+at)dudv 

—  oo  —  oo  \  r  t 


where 
A’ (u,v) 


°0  00 
—  00—00 


A  (u  V  ) 
l  /  l 


Sin  (u-u i  )  7T 
(u-Ui  )  TT 


SinJv-VjJjr^-i  (o-a  1 )  du 
(v-Vj  )  7T 


where  ai  =  a(u,v).  In  other  words  A'is  the  weighted 


average  of  values  of  A  over  neighbouring  wavelengths  and 


«; 
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directions.  This  is  true  of  a  square  region  S  with  side  2R, 
where  k  =  tt/R.  R  should  be  very  large  compared  with  the 
wavelengths  associated  with  most  energy  in  the  spectrum.  C, ' 
also  vanishes  outside  S.  Further  simplification  results  in 

P~P o  -  l/2gA  =  -Re  f  f  o2 A'A,e^^CJtdu  dv 

- CT  —oo  —  oo  — 

P  y 

Aj.  stands  for  A'  (-u,  -v) ,  and  is  the  amplitude  of  the  wave 
component  opposite  to  A  (u,v).  This  shows  that  1)  fluctu¬ 
ations  in  the  mean  pressure  p  arise  only  from  opposite  pairs 
of  wave  components  in  the  spectrum;  2)  that  the  contribution 
to  {p>  from  any  opposite  pair  of  wave  components  is  twice 
their  frequency  and  proportional  to  the  product  of  their 
amplitudes;  3)  and  that  the  total  pressure  fluctuation  is 
the  integrated  sum  of  the  contributions  from  all  opposite 
pairs  of  wave  components,  taken  separately. 

For  a  more  rigorous  treatment  of  the  subject,  the 
reader  is  referred  to  the  original  work  of  Longuet-Higgins . 

1.5  It  should  be  noted  here  that  the  earlier  work  on  micro¬ 
seisms  was  mainly  connected  with  the  oscillations  whose 
frequency  was  .25  to  .125  hz  or  lower.  It  was  not  until 
later  that  the  study  of  high  frequency  microseisms  was 
feasible  when  proper  instrumentation  was  made  available. 
Earlier  workers  were  also  contented  with  temporal  comparison 
of  microseisms.  An  attempt  was  made  first  in  Ottawa  (Macel- 
wane,  1952)  during  the  conferences  of  the  Eastern  Section  of 
the  Seismological  Society  of  America  to  correlate  and  map 
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microseisms  in  Canada  and  the  U.S.,  but  nothing  seems  to 
have  come  out  of  this  study,  at  .least  no  literature  can  be 
cited.  It  was  in  1952  that  a  programme  was  initiated 
by  the  Association  of  Seismology  of  the  International 
Geodetic  and  Geophysical  Union  to  record  simultaneously 
microseisms  of  the  whole  globe.  Like  its  predecessor  it 
did  not  seem  to  attract  the  attention  of  many  workers  or  at 
least  no  work  came  to  the  notice  of  the  present  writer. 

1.6  In  addition  to  the  temporal  and  spatial  correlations 
of  microseisms,  attempts  were  made  to  determine  the  direc¬ 
tion  of  approach  of  microseismic  waves  by  assuming  that 
microseisms  are  Rayleigh  waves.  Tromsdorff  (Macelwane,  1952) 
and  Ramirez  (1940)  independently,  using  the  tripartite  method, 
found  that  the  bearings  determined  from  arrival  time  dif¬ 
ferences  across  the  array,  pointed  towards  a  low  pressure 
area  over  the  sea.  Later  attempts  by  the  U.  S.  Naval  Aerolog- 
ical  Service  under  Gilmore  (1946)  in  tracking  cyclones  by  the 
Ramirez  method  did  not  prove  too  encouraging.  Jensen  (1958) , 
on  the  assumption  of  a  Rayleigh  wave  nature  of  microseisms, 
devided  a  vector  method  for  determination  of  direction  of 
approach.  He  mentions  that  at  the  moments  at  which  the  Z 
component  is  maximum,  the  slopes  of  the  two  horizontal 
records  are  read.  The  plot  between  these  two  slopes  gives  the 
velocity  of  the  earth  particle.  The  mean  of  several  velocity 
vectors  gives  the  direction  of  approach.  Jensen  (1958,  1959, 
1961,  1965)  using  this  method  determined  the  direction  of 
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approach  at  Copenhagen  and  at  several  places  in  Greenland. 

He  found  a  strong  correlation  between  the  direction  of 
approach  and  the  bearings  of  an  atmospheric  low  pressure 
area  over  the  North  Sea.  Darbyshire  (1954)  using  Lee's  method 
of  direction  of  approach  arrived  at  results  which  agreed 
with  meteorological  charts.  In  Lee's  method  (1935)  first  the 
seismograms  are  corrected  for  instrument  phase  response  and 
then  the  phases  are  read  off  the  record  for  each  component. 
Phase  differences  between  each  pair  of  the  three  components 
are  then  calculated.  These  phase  differences  are  used  to 
determine  the  nature  of  microseismic  waves.  Lee  (ibid)  con¬ 
cluded  that  all  microseisms  recorded  at  Kew  Observatory  were 
Rayleigh  waves.  On  the  other  hand  Darbyshire  (op.  cit.) 
showed  that  there  were  appreciable  love  waves  present  in  the 
microseisms  he  studied.  This  was  in  accord  with  observations 
made  by  Leet  (Lee,  1935).  In  Lee's  case  the  measurements  were 
very  few.  The  phase  response  for  each  instrument  was  cal¬ 
culated  at  a  few  discrete  periods  and  the  correction  was  made 
in  the  time  domain.  This  method  is  inferior  to  one  in  which 
the  overall  response  of  the  instrument  is  computed  and  then 
matched  to  the  measured  calibration  curve  of  the  instrument. 

In  none  of  the  papers  cited  there  is  ever  a  mention  of 
matching  of  the  observed  calibration  curves  against  calcu¬ 
lated  values  as  is  discussed  later  in  this  work.  A  more 
serious  source  of  error  lies  in  the  measurement  of  the  phases 
which  are  read  directly  from  the  traces.  Darbyshire  (op. cit.) 
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magnified  the  traces.  Even  with  magnification  inaccuracies 
are  bound  to  result  in  these  measurements. 

The  workers  on  the  North  American  continent,  how¬ 
ever,  did  not  find  the  determination  of  the  direction  of 
approach  so  simple.  Tripartite  stations  were  set  up  to  de¬ 
termine  the  direction  of  approach.  Considering  that  an 
equilateral  triangle  will  come  closest  to  giving  the  most 
accurate  results,  Dinger  (1952)  points  out  that  the  estimate 
of  direction  even  in  this  case  is  +'11°;  a  spread  of  22°.  The 
waves  often  come  from  different  directions  resulting  in  a 
very  low  coherency.  The  wave  trains  differ  not  only  in 
their  direction  of  arrival  but  in  period  and  wave  type. 

Thus,  as  will  be  confirmed  in  this  work,  a  simple  straight¬ 
forward  method  does  not  seem  to  yield  any  satisfactory  re¬ 
sults,  though  claims  contrary  to  this  do  exist  in  literature. 

1.7  The  most  intricate  problem  connected  with  microseisms 
seems  to  be  the  physical  nature  of  these  waves.  That  they 
are  not  purely  Rayleigh  waves  is  now  established.  In 
addition  to  the  fundamental  Rayleigh  mode,  the  existence  in 
microseisms  of  higher  Rayleigh  modes  and  also  of  different 
Love  modes,  has  been  reported  by  workers  such  as  Gutenberg 
(1958),  Ewing  et  al  (1957,  p  185  -  189).  Phinney  (1961),  in 
a  theoretical  consideration,  suggested  possible  existence  of 
leaking  P  modes  which  could  be  propagated  over  long  distances. 
Douze  (1964)  obtained  results,  using  seismometers  at  dif¬ 
ferent  depths,  which  point  out  that  noise  in  general  consists 
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of  an  arbitrary  combination  of  several  Rayleigh  modes. 

Several  authors  have  discussed  the  possible  existence  of 
body  waves  in  the  short  period  microseisms,  (Roden  (1965) , 
Seriff  et  al  (1965),  and  Gupta  (1965)).  Douze  (1967)  states 
"Previous  studies  of  short  period  seismic  noise  have  often 
assumed  that  only  surface  waves  are  present  in  short  period 
noise.  However,  the  experimental  results  obtained  to  date 
cannot,  in  general,  be  explained  in  this  fashion.  The  pres¬ 
ence  of  random  body  wave  noise  must  also  be  taken  into  con¬ 
sideration."  Having  examined  the  noise  in  the  range  6.0  - 
0.3  sec,  Douze  concludes  that  surface  waves  generally  pre¬ 
dominate  at  the  longer  periods  (of  the  period  range  discussed) , 
mostly  4-6  sec. ,  while  the  body  waves  appear  at  shorter 

i 

periods,  mostly  0.3-2  sec.,  at  quiet  sites.  He  also  states 
that  not  all  the  data  could  be  interpreted  to  define  the 
wave  type  present. 


CHAPTER  II 


METHOD  OF  ANALYSIS 

II. 1  The  aim  of  the  present  work  is  to  compare  the  micro- 
seismic  waves  at  different  stations  in  Canada  and  try  to 
identify  the  wave  types  present  in  the  data.  In  other 
words  more  emphasis  is  placed  on  the  mechanism  of  propagation 
of  microseisms  from  station  to  station  which  are  separated 
by  large  distances.  There  is  an  excellent  network  of  seismo¬ 
graph  stations  in  Canada  (see  Fig.  II. 1)  and  most  of  these 
have  the  facility  for  recording  on  three  short  period  and 
on  three  long  period  instruments.  A  casual  observation  of 
these  records  shows  existence  of  unmistakable  microseisms  on 
the  records.  Long  period  instruments  especially  show  the 
sinusoidal  type  of  motion  with  beat  phenomenon  superimposed 
on  it.  Winter  months  seem  to  be  the  time  when  microseismic 
activity  is  maximum  with  a  peak  activity  in  January  (Basham 
and  Whitham,  1966) . 

January  1965  records  showed  very  good  trains  of  these 
complex  waves  and  records  for  January  1-5  were  picked  out  at 
the  Leduc  observatory  some  20  miles  south  of  Edmonton.  Pre¬ 
liminary  work  showed  a  good  possibility  for  identification 
of  Rayleigh  and  Love  waves.  It  was  then  thought  worthwhile 
to  examine  the  records  at  some  of  the  other  operating  stations 
for  the  same  dates  and  times  and  to  attempt  a  temporal  and 
spatial  correlation  of  these  waves.  The  records  for  the 
same  time  were  obtained  from  Victoria,  Penticton,  Resolute 


' 

. 


Fig.  II.  1 


STANDARD  SEISMIC  STATIONS  IN  CANADA 
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and  Ottawa  stations.  All  of  them  showed  good  microseismic 
activity  for  the  period  mentioned.  The  next  step  then  was 
to  statistically  analyse  these  records  for  which  it  was 
necessary  to  digitize  the  seismograms.  The  seismogram  dig¬ 
itizer  as  developed  by  Wickens  and  Kallar  (1967)  at  the 
Dominion  Observatory,  Ottawa  was  used  for  this  work. 

II. 2  The  general  principle  of  the  digitizer  is  that  the 
seismogram  is  mounted  on  a  drum,  the  same  in  physical  di¬ 
mensions  as  the  recording  drums  used  in  the  network  stations, 
and  the  drum  is  driven  by  a  motor;  the  trace  is  followed 
optically  while  the  drum  is  rotating,  and; the  transverse 
movement  of  the  drum  required  to  follow  the  fluctuations  of 
the  trace  can  be  changed  into  equivalent  voltages  by  a  pre¬ 
cise  linear  potentiometer;  potentiometer  measurements  are 
then  fed  to  a  digital  voltmeter;  once  a  sample  is  held  in  the 
voltmeter  it  only  remains  to  punch  this  in  suitable  format 
onto  cards.  For  more  details  the  reader  is  referred  to  the 
paper  cited. 

For  the  purpose  of  this  work  the  voltage  sensitivity 
was  kept  at  .0325  volts/mm.  The  calibration  was  done  before 
each  day's  work  and  an  average  figure  was  taken. 

For  the  Edmonton  records  the  total  length  of  the 
record  digitized  was  five  minutes.  Edmonton  records  are  run 
twice  as  fast  as  those  at  other  stations.  The  total  length 
of  record  digitized  at  the  other  stations  was  twenty  minutes. 
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There  are  twelve  hundred  samples  in  each  case  so  that  the 
digitizing  interval  for  Edmonton  records  was  one  fourth  that 
of  other  stations,  the  actual  figures  being  0.25  sec  and 

I. 0  sec.  Thus  the  nyquist  frequency  for  Edmonton  records  is 

2.0  hz  and  that  for  the  records  at  other  stations  0.5  hz. 

The  shortest  period  estimate  in  these  two  cases  are  thus 

* 

0.5  sec  and  2  sec,  respectively. 

II.  3  Digitization,  or  sampling  of  data  at  discrete  inter¬ 
vals  places  a  restriction  on  the  data.  Statistical  estimates 
of  the  data  for  frequencies  above  1/2  At,  where  At  is  the 
sampling  rate,  do  not  exist  in  a  mathematical  form.  As  a 
result,  if  such  a  frequency  does  exist  in  the  data,  and 

if  we  call  it  (1/2  At  +  A) ,  then  it  will  be  folded  down  to 

(1/2 At  -  A),  since  Cos  (tt  +  6)  =  Cos  (g  -  6).  This  process 

is  called  aliasing.  The  frequency  (l/2At)  is  called  the 

"folding"  or  nyquist  frequency,  ^N.  The  frequencies  f, 

2f  -  f,  2f  +  f,  4f  -  f,  4f,  +  f,  and  so  on  are  called  the 
N  N  N  N 

"aliases"  of  one  another,  and  f  is  called  the  principal 

"alias".  The  "aliased  spectrum"  P  (f)  is  the  result  of 

a 

"aliasing"  P(f).  The  principal  part  of  the  aliased  spectrum 

P  (f)  is  the  part  of  P  (f)  which  corresponds  to  principal 
A  a. 

aliases,  positive  and  negative.  Thus  the  contributions  to  the 
power  pectrum  from  frequencies  greater  than  fN  cause  problems 
of  aliasing. 

Aliasing  is  an  unavoidable  consequence  of  digitiza- 
It  seems  desirable  to  evolve  a  workable  scheme  of 
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digitizing  the  data  in  some  definite,  but  not  at  equally 
distributed  intervals,  and  estimate  the  power  spectrum  with¬ 
out  aliasing.  No  such  scheme  seems  to  have  been  developed 
so  far. 

Aliasing  causes  more  serious  problems  with  data 
like  microseisms,  where  signal  to  noise  ratio  is  low.  Con¬ 
tributions  to  the  power  spectrum  from  high  frequency  noise 
present  in  microseisrns  because  of  aliasing  poses  an  important 
and  potential  threat  to  the  statistical  analysis  of  the  data. 
Identification  of  wave  types  may  become  difficult  because 
of  the  presence  of  high  frequency  noise  folded  into  the 
data  band. 

II. 4  If  one  is  looking  for  an  event  whose  frequency  is  much 
lower  than  the  folding  frequency,  one  can  safely  assume  that 
aliasing  problems  are  not  statistically  important  provided 
the  sampling  rate  is  sufficiently  high  so  that  negligible 
power  is  present  in  the  original  data  above  the  folding 
frequency.  Due  to  the  frequency  response  of  the  instru¬ 
ments  it  is  believed  that  this  is  not  a  serious  problem 
for  the  present  data.  As  mentioned  in  II. 2  the  lowest 
period  estimate  for  Edmonton  records  is  0.5  sec  and  the  same 
for  other  records  is  2.0  sec.  For  the  microseisrns  with  a 
peak  amplitudes  at  6  sec  and  4  sec,  the  statistical  estimate 
for  Edmonton  records  is  more  reliable  than  that  for  other 
stations  because  the  folded  power  will  be  less  serious  a 
problem  in  this  case.  An  example  of  the  original  trace 
from  the  seismogram  as  compared  with  that  of  the  digitized 
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trace  is  shown  in  figure  II. 2. 

II. 5  In  order  to  identify  the  wave  types  it  is  useful  to 
determine  the  cross  correlation  function  for  two  given  traces 
and  also  find  the  phase  difference  between  the  two.  Before 
making  this  calculation  it  is  necessary  to  determine  the 
phases  of  the  two  instruments  involved  and  make  the  necessary 
correction.  After  this  phase  correction  is  made,  any  cal¬ 
culation  of  phase  difference  between  the  two  instruments 
becomes  meaningful. 

There  are  two  ways  of  doing  this  phase  correction. 

In  the  first  method,  one  calculates  the  phases  from  the 
transfer  function  of  the  instrument  for  given  periods  and 
takes  the  phase  value  for  a  given  instrument  at  the  period 
one  is  interested  in.  In  the  case  of  microseisms  it  is 
around  6  sec.  Then  one  reads  the  phases  off  the  trace  and 
makes  the  necessary  correction  depending  on  whether  the  in¬ 
strument  is  leading  or  lagging.  The  disadvantages  in  this 
method  are: 

1.  The  trace  that  appears  on  the  seismogram  is  not 
•  exactly  a  "6  sec"  event.  It  is  composed  of 

various  other  frequencies.  Assuming  that  this 
trace  is  purely  of  one  given  frequency  is 
erroneous . 

2.  The  method  of  determining  the  phase  from  the 
trace  is  crude  and  inaccurate. 

In  the  second  method  one  calculates  the  transfer 
function  in  the  frequency,  domain.  Let  us  call  this  H  ( cu)  - 
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Then  one  takes  the  Fourier  transform  of  the  digitized  trace 
also  in  the  frequency  domain.  Let  us  call  this  F  (  go)  .  The 
division  of  F  ( co)  by  H(w)  in  the  frequency  domain  then  gives 
the  corrected  frequency  content  of  the  signal.  In  this 
method  phases  at  all  frequencies  are  taken  care  of.  The  in¬ 
verse  Fourier  transform  of  F(w)/H(w)  thus  puts  the  trace 
back  into  time  domain  with  phase  corrected.  In  the  series 
F(^),  the  first  term  is  the  d.c.  level  of  the  trace,  a 
meaningless  quantity  in  the  present  context.  The  next  two 
terms  correspond  to  periods  1024  sec  and  512  sec  for  Vic¬ 
toria,  Penticton,  Resolute  and  Ottawa,  and  256  sec  and  128 
sec  for  Edmonton.  The  effect  of  these  long  periods  is  re¬ 
moved  by  equating  the  Fourier  coefficients  for  these  frequen¬ 
cies  to  zero,  since  most  of  the  power  at  these  very  long 
periods  is  due  to  instrumental  drift. 


II. 6  It  is  now  necessary  to  examine  the  transfer  function 
for  the  galvanometer-seismometer  system  of  the  long  period 
instruments  used  in  recording  the  microseisms.  The  galvano¬ 
meter  is  a  low  pass  and  the  seismometer  is  a  high  pass  device 
in  terms  of  frequency.  The  velocity  sensitivity  for  the 
seismometer  can  be  easily  shown  to  be 
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and  the  velocity  sensitivity  for  the  galvanometer  is 

9  (s)  =  „  _ 1 _  (2) 

E(s)  *  1  +  2?g(s/wg)  +  (s/wG)2 


where  E(s)  is  the  voltage  developed  in  the  coil,  s=iw, ^ 

being  the  frequency,  co  and  m  the  natural  frequencies  of 

O  G 

the  seismometer  and  galvanometer,  r  and  r  are  the  damping 

o  G 

coefficients  for  seismometer  and  galvanometer,  0(s)  and  V(s) 
are  the  velocity  sensitivities  for  galvanometer  and  seis¬ 
mometer  respectively.  Multiplication  of  (1)  by  (2)  gives 

-/-?■  =  T_  (s)xt  (s)  ,  where  T0(s)  and  T_(s)  are  the  transfer 
V  iS)  b  b  b  b 

functions  of  seismometer  and  galvanometer  respectively  as 
given  by  (1)  and  (2) r  assuming  negligible  coupling  between 
the  two  systems. 

For  the  seismometer,  since  2tt/T,  from  (1) 

i)  as  T-*0,  E  (s)/V  (s)-*G,  a  constant  and  the  phase  difference 

is  0°.  Or 

ii)  as  hence  E(s)/V(s)+  -G(w/ms)2  a  line  of  slope 

2  on  a  log-log  plot  ie.  12db/octave,  and  the  phase 
difference  is  180°. 

Similar  results  hold  for  the  galvanometer  and  these  as- 
symptotic  values  are  shown  in  Figure  II. 3. 
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At  the  natural  frequency  of  the  seismometer  go  = 


GO 


hence 


E (s)  _  iG 
V (s )  2C* 


so  the  phase  is  90' 


E  ( s )  G 

i)  If  =  1, | vTsyl  =  2  ^own  6db  from  maximum  or 

ii)  using  -  the  common  condition  for  maximum  linear 

i  E  ( s )  G 

frequency  response,  then  |  -|s  j- 1  =  yj  down  3db  from  maximum 


II. 7  The  exact  determination  of  the  transfer  function  is 
not  possible  but  can  be  achieved  very  closely  by  comparing 
the  transfer  function  with  the  calibration  curves  measured 
routinely  at  the  various  observatories.  Calibrated  curves 
which  give  velocity  sensitivity  vs.  period  on  log-log  plot 


are  available  for  each  instrument  at  each  station.  A 
cursory  look  at  the  curves  indicate  that  there  is  no  ap¬ 
preciable  phase  difference  between  any  two  given  instru¬ 
ments  at  a  given  station.  A  typical  calibration  is  shown 
in  Figure  II. 4.  The  curve  in  this  figure  was  drawm 
approximately,  by  fitting  a  smooth  curve  through  the 
measured  points. 


The  transfer  function  as  given  in  II. 6  can  be 


written  as 


0  (s)  _  *  s 
E  (s )  a  a), 


1  +  2C 


S  GO. 


s_ 

go. 


1  +  2c 


G  u). 


s_ 

go" 


» 


station: 


OTTAWA 


<£«  45*23'38"N  X*  75M2'57"W  Altitude  83M 


Foundation.  Boulder  clay  on  limestone 


Dotes  of  Calibration  1 


COLUMBIA  l 

COLUMBIA  N-s  0  jrebt  7  _  i9g4 
COLUMBIA  C-W 


Fig.  II. 4 
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The  quantities  U)  ,  co  natural  frequencies  of  seis- 

O  G 

mometer  and  galvanometer  are  known.  s=i^o,  (co=27t/T)  can  be 

changed  at  regular  intervals  to  give  ^-|^yas  a  function  of 

s.  an<3  the  damping  coefficients  of  seismometer  and 

o  G 

galvanometer  can  be  varied  to  give  a  family  of  curves  of 

0(s)/E(s)  as  a  function  of  s.  For  each  and  £  there  will 

o  G 

be  one  such  curve.  We  have  the  measured  velocity  sensitivity 
at  a  series  of  frequencies  from  the  observatories.  These 
are  then  matched  with  the  family  of  curves  obtained  as  above. 
The  closest  matched  curve  as  illustrated  in  Figure  II. 5  is 
considered  as  the  velocity  sensitivity  response  to  period  and 
the  damping  coefficients  may  be  noted  for  this  curve.  In 
table  II. 1  are  enumerated  these  damping  coefficients  for 
all  stations  as  determined  from  matching  of  the  observed 
curve  with  the  calculated  curve.  Matching  does  not  neces¬ 
sarily  mean  matching  at  all  frequencies.  The  discussion  of 
the  transfer  function  for  seismometer  and  galvanometer  in¬ 
dicates  that  matching  has  to  be  done  only  in  the  region  be¬ 
tween  T^  and  T_  where  T_  and  T_  are  the  natural  periods  of 

S  VJ  OO 

the  seismometer  and  galvanometer  respectively.  Once  the 
portion  of  the  curve  between  the  two  points  is  closely 
matched,  the  gain  factor  A  is  computed  from  the  calibration 
curve.  Table  II. 2  gives  these  gain  factors.  Figure  II. 6 
shows  a  comparison  between  the  digitized  traces  and  the 
phase  corrected  traces.  Note  that  there  is  a  change  in 
amplitude  between  the  sets  of  traces  due  to  normalization 


» 

. 


Log  Velocity  Sensitivity 


RESPONSE  OF  THE  SEISMOMETER-GALVANOMETER  SYSTEM 


AT  OTTAWA 


FIG  II. 5 
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DAMPING  COEFFICIENTS 

TABLE  II. 1 


VERTICAL 

NORTH- 

•SOUTH 

EAST- 

-WEST 

Seis 

Galv 

Seis 

Galv 

Seis 

Galv 

Victoria 

1.4 

0.35 

1.4 

0.35 

1.4 

0.35 

Penticton 

1.4 

0.35 

1.4 

0.35 

1.4 

0.35 

Edmonton 

0.5 

0.7 

0.5 

0.7 

0.5 

0.7 

Resolute 

0.9 

.  35 

0.7 

0.7 

0.9 

0.35 

Ottawa 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

MAGNIFICATION  FACTORS 

TABLE  I I . 2 

)  . 

VERTICAL  NORTH- SOUTH 

EAST-WEST 

Victoria 

344 

480 

324 

Penticton 

275 

205 

180 

Edmonton 

390 

384 

372 

Ottawa 

1000 

1200 

1400 

Resolute 

3000 

1000 

3037 

DIGITIZED  TRRCE  Z 


DIGITIZED  TRRCE  N 


MVrW»WA'jt 


DIGITIZED  TRRCE  E 


0.00 


16.00 
H - 


qo 


.00 


TIME  IN  SECS  #  CX101  ) 


48.00 
— I - 


64.00  80.00  96.00 


l 


PHASE  C0RR.  TRACE  Z 


■■  Yw^yv'. 


PHASE  CORR.  TRRCE  N-S 

— ^V«^vvVwV^-Vsr*^V^o^r^\V>/^r^*# — 


PHASE  CORR.  TRRCE  E-W 


Fig.  II. 6 
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and  that  the  very  long  period  components  have  been  removed 
as  mentioned  previously. 

II. 9  The  information  available  for  the  present  study  lies 
in  the  digitized  results  for  the  three  components  Vertical, 
North-South  and  East-West  at  the  recording  sites  of  Victoria, 
Penticton,  Edmonton,  Resolute  and  Ottawa.  The  dates  on  which 
the  data  were  taken  are  January  1,  through  5  of  1965  when 
there  was  an  atmospheric  low  pressure  area  on  the  Pacific 
Ocean,  as  shown  in  figure  II. 7. 

Tables  II. 3-7  give  the  starting  time  of  digitization, 
the  number  of  "records"  digitized,  the  length  of  each  "record" 
and  the  interval  time  between  the  records.  From  these  tables 
it  can  be  seen  that  there  are  16  records  for  Victoria,  8  for 
Penticton,  39  for  Edmonton,  13  for  Resolute  and  8  for  Ottawa. 
Two  records  for  Victoria  and  one  for  Penticton  could  not  be 
analyzed  as  alpha-characters  were  punched  on  the  cards  in¬ 
stead  of  numbers  in  the  process  of  digitization.  All  the 
data  was  stored  on  tape  and  the  tape  was  subsequently  used 
for  the  whole  analysis. 

The  center  of  the  storm  did  not  remain  stationary 
during  the  time  January  1-5.  The  atmospheric  low  pressure 
area  on  the  Pacific  Ocean  changed  its  position  from  day  to 
day.  From  the  known  coordinates  of  the  storm  centre  and  the 
recording  sites,  distances  between  each  site  and  each  loca¬ 
tion  of  atmospheric  low  pressure  area  and  the  bearings  from 
one  to  the  other  were  calculated.  These  distances  and  bear- 
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ings  are  shown  in  Tables  II.8a-83. 

II. 8  Calculation  of  the  autopower  spectrum  of  a  stationary 
time  series  has  been  dealt  with  by  Weiner  (1964) ,  Blackman  and 
Tukey  (1958) ,  and  Lee  (1964)  among  others.  The  methods  have 
been  standardized  and  used  by  geophysics  researchers  as  rou¬ 
tine  procedures.  A  discussion  of  these  methods  is  avoided  in 
this  work  because  of  the  availability  of  more  exhausitve 
treatments  elsewhere.  Suffice  it  to  say  that  Fourier  analysis 
technique  is  exploited  to  obtain  the  power  of  the  series.  The 
Fourier  transform  of  the  series  is  first  calculated,  a  filter¬ 
ing  or  smoothing  if  necessary  is'  done  in  the  frequency  domain, 
correlation  techniques  are  performed  in  the  lag  domain  and 
finally  inverse  Fourier  transform  of  the  series  gives  the 
power  as  calculated  at  discrete  values  of  period  or  fre¬ 
quency.  The  basic  assumption  is  that  stationary  time  series 
is  composed  of  a  large  number  of  sinusoidal  waves  varying 
in  frequency.  The  Fourier  analysis  of  the  data  calculates 
the  power  at  all  those  frequencies  present  in  the  data  which 
are  below  the  nyquist  frequency.  The  peaks  in  the  power 
spectrum  correspond  to  the  "signal"  otherwise  obscured  by 
noise  in  the  time  series.  The  same  techniques  can  be  used 
for  two  series  to  obtain  the  power  between  them.  In  one 
case  the  power  is  called  "autopower"  of  the  series  and  in 
the  second  case  it  is  called  "crosspower"  between  the  two 
series.  Crosspower  calculation  also  has  an  option  whereby 
phases  and  coherencies  between  two  given  series  are  cal- 
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culated.  Crosspower  is  more  effectively  used  in  detecting 
a  signal  common  to  both  series.  Thus  it  is  helpful  in  de¬ 
tecting  any  kind  of  phenomenon  which  is  recorded  on  two 
instruments  simultaneously. 

II. 10  Microseisms  can  be  assumed  to  be  stationary  in  time 
if  the  record  is  long  compared  with  the  period  of  the  event 
concerned.  Inherent  in  the  assumption  is  the  point  that  it 
is  possible  to  predict  the  statistical  behaviour  of  the  time 
series  at  a  time  T,  if  its  behaviour  is  known  at  some 
earlier  time.  The  general  procedure  adopted  by  the  workers 
in  the  field  has  been  directed 

1)  to  identify  the  wave  type. 

2)  to  determine  the  direction  of  approach  of 
microseisms . 

3)  to  find  the  effect  of  the  geological  features 
on  their  propagation. 

4)  to  investigate  the  possibility  of  tracking  down 
the  storms  on  the  ocean. 

The  present  work  deals  mainly  with  the  mechanism 
of  propagation  of  microseisms  over  a  continent  and  an  attempt 
at  the  identification  of  different  wave  types  present  in 
microseisms.  An  attempt  is  also  made  to  determine  the  direc¬ 
tion  of  approach.  Since  the  data  is  available  on  three  com¬ 
ponents  at  each  station  at  different  times,  the  crosspower 
analysis  between  vertical  and  horizontal  components  would 
yield  information  on  the  presence  or  otherwise  of  Rayleigh 
wave,  and  that  between  two  horizontals  would  give  the  same 
information  about  Love  waves. 
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CHAPTER  III 


POWER  SPECTRA  OF  MICROSEISMS 

III. I  Autopower  spectrum  was  calculated  for  a  total  number 
of  seventy-nine  records,  thirteen  of  which  were  taken  from 
Victoria,  seven  from  Penticton,  thirty-nine  from  Edmonton, 
eight  from  Ottawa  and  twelve  from  Resolute.  In  each  record 
there  are  three  traces,  viz.  one  vertical  and  two  horizontal 
components.  Three  sets  of  crosspower  were  also  calculated 
along  with  the  coherencies  and  phases,  one  between  vertical 
and  east-west  components,  a  second  between  vertical  and  north- 
south  components  and  a  third  between  the  two  horizontal  com¬ 
ponents,  east-west  and  north-south.  Figure  III.l  is  a  sample 
of  autopower  spectrum  as  a  function  of  period.  In  Fig.  III. la 
the  autopower  spectrum  for  Edmonton  is  shown  and  in  Fig.  III. lb 
the  same  is  shown  for  Victoria. 

At  Edmonton  three  peaks  on  the  vertical  component  at 
2.3  sec,  3.3  sec  and  8.0  sec  and  two  on  the  east-west  and  north- 
south  components  at  2.5  sec  and  6.4  sec  can  be  observed.  The 
long  period  event  does  not  peak  either  at  8.0  sec  or  at  6.4 
sec  on  all  records;  the  period  varies  from  record  to  record, 
the  range  of  variation  being  5. 5-8.0  sec.  Similarly  the  3.3 
sec  event  varies  in  time,  variation  being  3. 3-4.9  sec.  On 
some  records  these  two  peaks  vanish  but  a  new  peak  at  5.2  sec 
is  clearly  shown,  The  short  period  event  is  always  present. 
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LOG  AUTO-POWER 


AUTO-POWER  SPECTRUM  OF  MICROSEISMS  AT  EDMONTON 
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its  range  of  variation  in  time  being  2. 3-2. 6  seconds.  The 
power  in  the  short-period  event  is  equal  to  or  greater  than 
the  power  in  the  longer  period  events. 

At  Victoria  the  6.4  sec  peak  is  present  on  all  three 
components.  A  peak  around  2.2  sec  can  be  seen  on  the  two 
horizontal  components  along  with  another  at  4.0  sec  for  the 
example  shown  which  also  appears  on  the  vertical  component. 

The  period  varies  in  time,  the  range  of  variation  being  5.4- 
6.6  sec,  4. 0-4. 6  sec  and  2. 0-2. 2  sec  for  the  three  events 
mentioned.  Unlike  at  Edmonton,  more  power  is  concentrated  in 
the  6.0  sec  event  than  in  the  other  two.  The  power  in  the 
shortest  period  event  is  almost  negligible  compared  to  the 
power  in  the  two  long  period  events  on  most  records.  The  spec¬ 
tral  behavior  of  microseisms  at  Penticton,  Resolute  and  Ottawa 
is  very  similar  to  that  at  Victoria. 

This  marked  contrast  between  Edmonton  and  the  other 
stations  in  the  behavior  of  spectra  especially  with  respect 
to  the  2.3  sec  event  presumably  is  related  to  the  contrast  in 
crustal  layering  at  these  recording  sites.  The  thick  sedim¬ 
entary  section  at  Edmonton  probably  accentuates  the  short- 
periods  . 

III. 2  In  figures  III. 2-8  is  shown  the  temporal  variation  of 
autopower  spectra  for  the  three  components  at  each  recording 
station.  The  upper  curve  in  these  diagrams  shows  the  vari¬ 
ation  of  magnitude  of  power  in  time,  the  lower  curve  shows 
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AUTO-POWER  PEAKS  AT  PENTICTON 
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AUTO-POWER  PEAKS  AT  EDMONTON  NORTH-SOUTH  COMPONENT 


m 


o 


H3M0d  DOT 


o 

i 


SQN003S  NI  aOIH3d 


FIG  III 


AUTO-POWER  PEAKS  AT  EDMONTON  EAST-WEST  COMPONENT 
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variation  of  period  in  time.  The  table  III.l  shows  the 
minimum  and  maximum  values  of  period  at  different  stations 
and  their  geometric  mean.  Not  much  statistical  importance 
can  be  attached  to  this  mean  value,  but  it  is  useful  as  a 
preliminary  guide  line  for  identifying  the  event.  If  we  now 
take  the  following  relationship  between  the  period  and  veloc¬ 
ity  viz.  T  =  0.4V  (Defant,  1961),  V  being  the  velocity  of  the 
wind  in  the  ocean  storm  center  we  can  easily  calculate  the 
high  and  low  values  this  velocity  attains  and  a  mean  value  of 
this  can  be  taken.  The  mean  velocity  as  calculated  from  the 
known  value  of  the  period  varies  between  15.5  and  17  knots. 

As  can  be  seen  from  meteorological  charts,  this  is  in  good 
agreement  with  the  value  of  wind  velocity. 

Table  III. 2  now  lists  the  periods  at  which  there  exists 
an  autopower  peak  for  each  component  and  also  the  periods  at 
which  there  exists  a  crosspower  peak  between  two  given  com¬ 
ponents  at  each  station.  The  values  in  this  table  are  actually 
the  mean  values  calculated  before  as  shown  in  table  III.l. 

III. 3  Some  salient  features  of  the  histograms  and  curves 
shown  in  figures  III. 2-8  can  be  cited  as  follows: 

1.  There  are  random  changes  in  the  period  with 
time. 

2.  At  a  given  time  and  a  given  site  the  period 
at  which  peak  is  observed  in  the  power  spec¬ 
trum  of  one  component  does  not  always  corres¬ 
pond  to  the  period  at  which  a  peak  is  observed 
on  another  component.  This  discrepancy  varies 
from  a  fraction  of  a  second  to  as  much  as  3 
seconds. 
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TABLE  III.l 

AUTOPOWER  AND  CROSS  POWER  PEAK 
RANGES  IN  TERMS  OF  PERIOD  IN  SECONDS 
AND  WIND  VELOCITIES  IN  KNOTS 
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TABLE  III. 2 


MICRO SEISMIC  PERIODS  AT  WHICH 


THERE  IS 

PEAK 

IN  THE 

POWER 

SPECTRUM 

Component 

Z 

N 

E 

ZN 

NE 

EZ 

Station 

Victoria 

5.9 

6.1 

5.8 

6.2 

6.4 

6.4 

Penticton 

6.2 

6.2 

6.2 

6.3 

6.2 

6.3 

Resolute 

6.3 

6.4 

6. 5 

6.4 

6.3 

6.5 

Edmonton 

6.8 

6.8 

6.2 

6.0 

6.8 

6.8 

Ottawa 

* 
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6 . 6 

6.1 
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r  , 
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3.  The  periods  can  not  be  correlated  spatially 
to  better  than  about  one  or  two  seconds. 

4.  At  some  stations  only,  there  is  a  rough 
correlation  between  changes  in  the  period  of 
maximum  microseism  power  and  the  changes  in 
the  wind  velocity  in  the  ocean  storm. 

5.  The  variations  in  period  are  remarkably  large 
at  Edmonton  compared  to  the  other  stations. 

The  variation  of  magnitude  of  power  is  almost 
oscillatory  in  time.  No  straightforward 
explanation  is  available  for  this  phenomenon. 

6.  The  power  reaches  a  maximum  value  and  then 
decreases  at  Victoria,  Penticton  and  Resolute. 

7.  At  Edmonton  no  such  systematic  variation  is 
evident. 

8.  At  Ottawa  the  variation  in  period  with  time 
and  the  variation  in  the  magnitude  of  the 
autopower  with  time  have  a  general  trend  of 
increase  throughout  the  time  interval  sampled 
by  the  data  sets,  suggesting  that  stable  storm 
conditions  have  not  yet  been  reached.  This 
may  be  due  to  another  storm  approaching  the 
East  Coast  at  the  same  time  the  storm  under 
consideration  is  receding  away  from  the  West 
Coast. 


III. 4  At  all  recording  sites  except  at  Edmonton  the  power 
reaches  a  maximum  value  on  all  three  components  at  approxim¬ 
ately  the  same  time.  The  behavior  of  power  as  a  function  of 
time  is  very  similar  for  all  three  components  at  these  four 
stations.  This  shows  that  the  mechanism  of  propagation  of 
microseisms  is  in  general  very  similar  in  the  portions  of  the 
crust  separated  by  these  four  sites.  If  the  microseisms  con¬ 
tain  both  Rayleigh  and  Love  waves,  which  seems  probable  from 
the  power  recorded  on  the  three  components,  then  both  these 
wave  types  seem  to  arrive  simultaneously  at  these  four  sites. 
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VARIATION  OF  AUTO-POWER  WITH  DISTANCE 
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This  feature  is  borne  out  by  the  statement  made  in  the  first 
sentence  of  this  section.  No  such  conclusion  can  be  reached 
on  Edmonton  data  because  there  is  a  marked  irregularity  at 
this  location  from  component  to  component  in  the  variation 
of  power  as  function  of  time. 

There  is  more  power  in  the  vertical  component  at 
Victoria,  Penticton  and  Ottawa  than  is  present  in  the  two 
horizontal  components.  At  Resolute  the  north-south  component 
registers  more  power  than  the  vertical  and  the  east-west  com¬ 
ponents.  This  presumably  has  something  to  do  with  the  crustal 
structure  at  these  sites.  For  a  Rayleigh  wave  travelling 
in  a  homogeneous,  isotropic  half-space  the  amplitude  ratio 

between  the  vertical  component  and  the  radial  component 

) 

should  be  1.5  and  hence  the  power  ratio  between  these  two  com¬ 
ponents  should  be  2.25.  Table  III. 3  shows  the  source  is 
neither  east-west  nor  north-south  for  any  recording  site  if  we 
assume  this  particular  power  ratio  to  be  true.  The  situation 
is  further  complicated  by  the  fact  that  Rayleigh  wave  motion 
is  not  the  same  for  a  multi-layered  medium  as  it  is  for  a 
half-space,  and  the  crustal  layering  at  the  various  sites  is 

markedly  different,  and  by  the  probable  presence  of  body  waves 
in  the  data. 

III. 5  Figure  III. 9  shows  the  spatial  variation  of  the  peak 
value  of  the  autopower.  There  is  a  maximum  value  of  power  at 
Penticton  but  a  decrease  on  either  side  of  this  station.  Why 
the  power  at  Victoria  is  smaller  than  that  at  Penticton  can 


_ 
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not  be  answered  at  all  if  we  assume  that  main  microseism 

signal  has  an  amplitude  which  decreases  as  —  and  that  the 

r11 

source  is  to  the  west  of  these  two  locations.  The  value  of 
n,  however,  can  be  calculated  from  the  power  recorded  at 
Penticton,  Edmonton  and  Ottawa  sites.  Mean  value  of  n  so 
calculated  is  2.1  for  the  vertical  component,  2.6  for  the 
north-south  component  and  also  2.6  for  the  east-west  compon¬ 
ent.  The  autopower  on  the  vertical  component  thus  varies  as 
— „  ,  and  hence  the  amplitude  varies  as  —  approximately,  as 

•  -L  JC 

r 

would  be  expected  for  a  simple  surface  wave.  If  one  were  to 

take  the  actual  figures,  one  would  observe  that  the  amplitude 

varies  as  —  where  n  is  between  1  (for  surface  waves)  and  2 
r 

(for  body  waves) ,  indicating  some  kind  of  a  quasi-surface  wave 
(perhaps  a  leaking  mode) . 

III. 6  The  velocity  sensitivity  spectra  for  three  stations 
Victoria,  Penticton  and  Edmonton  are  shown  in  Figures  III. 10, 
11  and  12.  At  all  these  sites  these  spectra  show  sequential 
peaks. 

Andrianova  et  al.  (1967)  in  a  treatise  published  in 
Seismic  Love  waves  (p.45)  point  out  that  for  foci  located  at 
the  depth  of  30-40  km  it  is  possible  to  record  in  a  crust 
of  30-40  km  not  only  the  first,  but  the  second  and  third  modes 
of  Love  waves  for  periods  less  than  7-8  sec.  The  dispersion 
curves  for  second  and  third  modes  in  a  crust  with  a  sedimen¬ 
tary  layer  are  marked  by  maxima  and  minima,  (p.  39,  Fig.  .15b, 
Model  8,  Andrianova,  op.  cit.)  According  to  the  same  authors 
specific  vibrations  with  quasi— stationary  periods  (less  than 
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VELOCITY  SENSITIVITY  SPECTRA  AT  VICTORIA 


i _ i _ i _ i _ i 

o 


00 

o 


_l _ J _ J _ I- 

O  vO 


J  I  L 


00 


I _ J _ 1_ 


OH  O  OH 

AIi  IA  I  iL  I  S  N3  S  AIiIOOTSA 


CD 

o 


o 

o 


DOT 


o 

i — I 


H 

H 

H 

CD 

H 

Ph 


VELOCITY  SENSITIVITY  SPECTRA  AT  PENTICTON  AND  EDMONTON 


vd  o  o  r-~  ^ 

.  •  •  •  • 


(N 


O 


O 


CO 


CM 


CO 

o 


o 


CO 

CO 

• 

Q 

1 — 1 

2 

O 

U 

w 

CO 

CM 

2 

• 

H 

rH 

Q 

O 

H 

Pi 

CO 

w 

• 

cu 

o 

e> 

o 

o 


CO 


CD 

o 

o 


AJ.IAII.I  SN3S  AiLI0013A  901 


FIGS  III. 11-12 


41 


2  secs)  and  long  duration  can  be  generated  corresponding  to 
resonance  in  the  system  consisting  of  a  surface  layer  and  a  • 
wave  guide  in  the  crust.  If  we  now  assume  that  the  Longuet— 
Higgins  theory  is  true  to  the  extent  that  the  ocean  water 
would  excite  elastic  waves  in  the  ocean  bottom  and  that  the 
oceanic  crust  acts  as  a  wave. guide  for  the  seismic  waves  to 
be  propagated  to  the  continent  and  thence  along  a  higher 
velocity  wave  guide  beneath  the  granitic  layer  then  we  satis¬ 
fy  the  first  condition  in  Andrianova's  theory,  that  the  source 
is  at  a  depth  of  around  30-40  km,  by  the  time  we  observe  the 
waves  on  the  continent.  Observe  that  the  crustal  thickness 
at  Victoria  is  in  excess  of  40  km. 

Thus  in  terms  of  Andrianova's  theory  Victoria  data 
could  be  interpreted  as  consisting  of  higher  modes  of  Love 
and  perhaps  of  Rayleigh  Waves,  since  the  vertical  component 
also  shows  sequential  peaks  in  the  velocity  sensitivity 
spectra.  It  would  be  however  impossible  to  interpret  Pen¬ 
ticton  and  Edmonton  data  in  terms  of  higher  modes  of  surface 
waves  as  these  modes  would  attenuate  at  larger  distances.  In 
addition  to  this,  a  sedimentary  layer  of  thickness  5  kms  with 
a  transverse  wave  velocity  of  2.6  km/sec  is  required  for  the 
2  sec  microseisms  to  be  recorded  at  Victoria.  In  view  of  the 
fact  that  such  layer  is  absent  at  Victoria,  and  in  view  of 
the  difficulty  of  interpreting  Penticton  and  Edmonton  data 
in  terms  of  higher  modes  of  surface  waves,  one  would  have  to 
look  for  a  different  interpretation  of  these  spectra. 


, 

. 


' , 
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Phinney  (1964)  working  on  earthquake  body  waves  of 
periods  5-50  secs,  obtained  curves  for  spectral  ratios  vs 
frequency  in  which  "characteristically  the  spectral  response, 

T  /  shows  sequence  of  peaks  due  to  resonance  effects  in  the 
crust."  Phinney  (1961)  also  mentions  in  a  theoretical  con¬ 
sideration  on  leaking  modes  "...  pressure  variations, 

.  .  .  due  to  standing  water  waves,  will  excite  this  mode 

(late  arriving  quasi-standing  mode)  at  least  as  efficiently 
as  they  excite  undamped  Rayleigh  waves.  The  noise  may  then 
leave  the  source  region  either  by  horizontal  propagation  in 
the  layer  or  by  means  of  the  P  waves  leaking  into  the  bottom. 
The  former  must  be  regarded  as  short-range  mechanism  where  as 

the  latter  would  carry  the  energy  to  appreciable  distances  as 

> 

body  waves."  We  merely  mention  these  two  statements  from 
Phinney  to  suggest  that  body  waves,  perhaps  leaky  modes,  are 
present  in  the  data. 

III. 7  Figures  III. 13-19  show  the  crosspower  spectral  abun¬ 
dances  at  each  station.  These  histograms  show  exactly  the 
same  features  as  those  mentioned  under  the  autopower  spectra. 
The  peak  shifts  in  period  from  record  to  record  at  a  given 
station,  the  shift  being  as  small  as  a  fraction  of  a  second 
and  as  large  as  3  seconds.  The  following  features  can  be  re¬ 
cognized  from  these  histogram: 

1.  On  Victoria  records  the  range  of  variation 

for  the  period  at  which  the  crosspower  peaks 
is  5. 6-6.9  seconds  for  the  vertical  and  North- 
South  components,  6. 1-6. 7  seconds  for  the 
North-South  and  East-West  components  and  5.9- 
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CROSS-POWER  PEAKS  AT  PENTICTON 
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CROSS-POWER  PEAKS  BETWEEN  VERTICAL  AND  NORTH-SOUTH  COMPONENTS  AT  EDMONTON 
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FIG  III. 15 


CROSS-POWER  PEAKS  BETWEEN  NORTH-SOUTH  AND  EAST-WEST  COMPONENTS  AT  EDMONTON 
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FIG  III. 16 


CROSS-POWER  PEAKS  BETWEEN  EAST-WEST  AND  VERTICAL  COMPONENTS  AT  EDMONTON 
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CROSS-POWER  PEAKS  AT  RESOLUTE 
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CROSS-POWER  PEAKS  AT  OTTAWA 
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FIG  III. 19 
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TABLE  III. 3 

RATIOS  OF  VERTICAL  TO  HORIZONTAL  POWERS 


Station  Z 

Power  EW 

Power 

Ratio 

NS  Power 

Ratio 

Victoria 

20 

15 

1.3 

15 

1.3 

Pentiction 

63 

25 

2.5 

40 

1.5 

Edmonton 

40 

16 

2.5 

9 

4.5 

Resolute 

0.4 

0.18 

2.2 

2.5 

0.16 

Ottawa 

2.5 

0.7 

3.5 

0.22 

10 

TABLE  III. 4 


POWER  VS  DISTANCE  (POWER  1_  ) 

rn 


Station 

Distance 

Power 

in  km 

Z 

N 

Penticton 

1111 

66 

40 

Edmonton 

1534 

40 

9 

Ottawa 

4344 

2.5 

.  84 

• 

Mean 

E 

n  for  Z 

n  for  NS 

n  for 

25 

1.6 

2.1 

1.3 

16 

2.7 

2.3 

3.4 

.23 

2.2 

4.4 

4.1 

2.1 

2.6 

2.6 
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7.3  seconds  for  the  east-west  and  vertical 
components . 

2.  On  Penticton  records  the  respective  figures 
are  6.1-6. 6  seconds,  5. 7-6. 6  seconds  and 

6 . 1-6 . 7  seconds . 

3.  On  Edmonton  records  the  peaks  are  more  scat¬ 
tered  in  period,  the  variation  being  5-8 
seconds,  a  spread  of  three  seconds  for  the 
long  period  event. 

4.  On  Resolute  records  the  corresponding  figures 
are  6. 0-6.9  seconds,  5. 8-6.9  seconds  and  6.2- 
6.9  seconds. 

5.  On  Ottawa  records  they  are  5. 6-7.1  seconds, 

5. 8-7.1  seconds,  and  5. 9-7.1  seconds  respec¬ 
tively. 

The  crosspower  spectra  thus  do  not  conclusively  point 
toward  a  well  defined  wave  type  occurring  within  a  range  of 

fraction  of  second  in  period.  The  data  is  scattered  from 

) 

trace  to  trace  at  almost  all  stations,  and  there  is  a  scatter¬ 
ing  from  station  to  station.  This  is  not  due  to  the  recording 
but  due  to  the  arriving  micro seisms  themselves  as  can  be  seen 
from  the  following  remarks  and  the  later  discussion  of  syn¬ 
thetic  microseisms  in  chapter.  IV.  The  coherency  values  rarely 
exceed  0.7  at  these  peak  values,  especially  where  the  vertical 
component  is  involved.  The  phase  between  the  vertical  and  a 
horizontal  component  is  rarely  -90  as  would  be  expected  of  a 
Rayleigh  wave  recorded  on  these  two  components.  The  phase 
between  the  two  horizontals  is  also  rarely  0°  or  180°  as  would 
be  expected  of  a  Love  wave  that  would  be  recorded  on  the  two 
horizontal  components.  Fig.  III. 20  shows  a  typical  example 
for  an  Edmonton  trace  of  a  plot  of  phase  vs  period  and  coher- 
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FIG  III. 20 
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ency  vs  period.  Phase  changes  rapidly  from  one  frequency 
to  the  other.  In  addition,  the  peaks  on  crosspower  vs  period, 
and  coherency  vs  period  are  not  at  all  correlated.  Not  one 
trace  at  all  these  stations  could  be  found  which  would  show 
all  the  three  quantities  correlated.  Coherencies  are  low 
because  of  the  lack  of  correlation  on  the  autopower  peaks  for 
various  components.  In  general  the  phase  also  varies  in  a 
random  way  and  seldom  takes  on  consistent  values  related  to 
any  simple  wave  type.  There  are  only  a  few  examples  of  well 
defined  Rayleigh  and  Love  wave  recordings  at  these  stations 
and  they  will  be  discussed  in  the  next  chapter. 


CHAPTER  IV 


WAVE  TYPES  IN  MICROSEISMS 

IV. 1  In  tables  IV. 1-5  are  listed  some  of  the  more  im¬ 
portant  events  recorded.  In  these  tables  are  shown  the 
period  at  which  either  there  is  a  peak  in  coherency  or  in 
crosspower,  and  the  values  of  the  coherency,  phase  and  the 
components  involved.  For  a  Rayleigh  wave  coming  from  the 
west,  the  coherency  between  the  east-west  and  vertical  com¬ 
ponents  has  to  be  high,  the  phase  between  the  two  has  to  be 
-90°  for  the  simplest  wave  type,  and  there  should  be  a  clear 
cut  peak  on  the  crosspower  spectra  between  these  components. 
For  a  Rayleigh  wave  arriving  from  the  north  the  two  components 
involved  are  vertical  and  north-south.  In  the  last  column 
of  these  tables  an  attempt  is  made  to  identify  the  wave  type 
on  this  basis. 

It  is  quite  obvious  that  many  spectral  peaks  could 
not  be  readily  identified  as  to  wave  type.  Some  of  these 
events  were  identified  as  Rayleigh  or  Love  or  the  first  mode 
of  Rayleigh  on  the  basis  of  phase  and/or  coherency  and/or 
peak  on  the  crosspower  spectrum.  An  event  with  a  peak  on 
coherency  and  crosspower,  and  a  phase  of  -90  y20  between 

the  vertical  and  a  horizontal  component  could  be  identified 
as  Rayleigh.  In  this  case  when  occasionally  the  phase  is 
approximately  —90^,  the  other  two  criteria  are  seldom  sat— 


. 


TABLE  IV. 1 
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VICTORIA 


Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh. 

ss  -power 

3 

5 . 8 

0.45 

-72 

EZ 

Yes 

Yes 

Rayleigh 

4 

6.7 

0.18 

-81 

EZ 

No 

Yes 

5.1 

0.52 

148 

EZ 

Yes 

Yes 

4.4 

0.47 

-73 

EZ 

Yes 

No 

5 

6.1 

0.41 

18.0 

NE 

Yes 

Yes 

Love 

6.6 

0.29 

107 

EZ 

No 

Yes 

Rayleigh 
(First  Mode) 

4.1 

0.43 

-3.9 

NE 

Yes 

Yes 

Love 

6 

6.4 

0.16 

104 

EZ 

Yes 

Yes 

Rayleigh 
(First  Mode) 

4.5 

0.47 

-103 

ZN 

Yes 

Yes 

Rayleigh 
(First  Mode) 

6.7 

0.55 

-49 

NE 

Yes 

Yes 

Love  (?) 

7 

6.0 

0.36 

-102 

ZN 

Yes 

Yes 

Rayleigh 

5.8 

0.24 

-70 

NE 

Yes 

Yes 

6.1 

0.22 

104 

EZ 

Yes 

Yes 

4.2 

0. 34 

-120 

ZN 

Yes 

Yes 

8 

6.0 

0.16 

-92 

EZ 

Yes 

Yes 

Rayleigh 

4.5 

0. 39 

180 

ZN 

Yes 

Yes 

6.1 

0.12 

8.2 

NE 

Yes 

Yes 

Love  (?) 

5.6 

0.24 

135 

NE 

Yes 

Yes 

Love  (?) 

9 

6.6 

0.40 

-131 

EZ 

Yes 

Yes 

6.0 

0.51 

7.0 

EZ 

Yes 

Yes 

5.7 

0.49 

-10 

ZN 

Yes 

Yes 

6.7 

0.47 

17 

NE 

Yes 

Yes 

Love  (?) 

c 

. 

■ 
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Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh . 

ss-  power 

10 

6.1 

0.37 

11 

EZ 

Yes 

Yes 

5.2 

0.62 

82 

EZ 

Yes 

Yes 

6.2 

0.28 

-17 

NE 

Yes 

Yes 

11 

6.1 

0.34 

1.8 

EZ 

Yes 

Yes 

6.1 

0.35 

100 

NE 

Yes 

Yes 

5.7 

0.59 

-0.9 

ZN 

Yes 

Yes 

12 

5.8 

0.23 

8.7 

EZ 

Yes 

Yes 

6.6 

0.35 

155 

NE 

Yes 

Yes 

Love 

6.6 

0.47 

9 

EZ 

Yes 

Yes 

4.9 

0.41 

-111 

EZ 

Yes 

Yes 

/ 

Rayleigh 

13 

6.1 

0.78 

4.4 

NE 

Yes 

Yes 

Love 

6.1 

.88 

-10 

EZ 

Yes 

Yes 

4.5 

.92 

4 

NE 

Yes 

Yes 

Love 

4.6 

.94 

-8.6 

EZ 

Yes 

Yes 

4.6 

.90 

-3.5 

ZN 

Yes 

Yes 

14 

6.1 

0.80 

-1.3 

NE 

Yes 

Yes 

Love 

6.1 

0. 86 

-0.3 

EZ 

Yes 

Yes 

6.2 

0.90 

1.0 

ZN 

Yes 

Yes 

4.6 

0.93 

12.0 

ZN 

Yes 

Yes 

4.6 

0.94 

4.0 

NE 

Yes 

Yes 

Love 

4.6 

0.93 

-5.1 

EZ 

Yes 

Yes 

15 

6.4 

0.73 

167 

NE 

Yes 

Yes 

Love 

6.2 

0.75 

-28 

EZ 

Yes 

Yes 

6.1 

0.75 

20 

4.6 

0.85 

20.5 

ZN 

Yes 

Yes 

4.6 

0.91 

-5.7 

NE 

Yes 

Yes 

Love 

4.6 

0.95 

-17 

EZ 

Yes 

Yes 

V 

, 

■ 

Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  cro¬ 

iod 

ency 

onent 

Coh. 

ss-power 

16 

6.4 

0.69 

-32 

EZ 

Yes 

Yes 

5.  8 

0.81 

2.2 

NE 

Yes 

Yes 

Type 


Love 


TABLE  IV. 2  50 

PENTICTON 


Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh. 

ss  -power 

18 

6.7 

0.45 

-76 

EZ 

Yes 

Yes 

Rayleigh 

6.7 

0.44 

10 

NE 

No 

Yes 

Love 

4.4 

0.83 

-1.7 

NE 

Yes 

Yes 

Love 

19 

6 . 6 

0.62 

9.9 

EZ 

Yes 

No 

7.1 

0.50 

-5.5 

NE 

Yes 

No 

Love 

6.1 

0.41 

6.3 

NE 

Yes 

Yes 

Love 

4.6 

0.88 

-5.8 

NE 

Yes 

Yes 

Love 

6.6 

0.65 

-13 

EZ 

Yes 

Yes 

20 

6.2 

0.48 

-111 

EZ 

Yes 

Yes 

Rayleigh 

6.2 

0.38 

-20 

NE 

Yes 

Yes 

4.6 

0.85 

0 

NE 

Yes 

Yes 

Love 

21 

6.1 

0.65 

15 

EZ 

Yes 

Yes 

6.4 

0.50 

-43 

NE 

Yes 

Yes 

4.5 

0. 85 

5 

NE 

Yes 

Yes 

Love 

22 

6.1 

0.21 

-50 

EZ 

Yes 

Yes 

6.4 

0.55 

-6 

NE 

Yes 

Yes 

Love 

4.6 

0.85 

11 

NE 

Yes 

Yes 

Love 

4.7 

0.51 

-76 

EZ 

No 

No 

Rayleigh 

23 

6.4 

0.64 

-3 

NE 

Yes 

Yes 

Love 

4.6 

0.91 

0.1 

NE 

Yes 

Yes 

Love 

4.9 

0.62 

-82 

EZ 

Yes 

Yes 

Rayleigh 

24 

6.0 

0.68 

-1 

NE 

Yes 

Yes 

Love 

4.5 

0.92 

-3 

NE 

Yes 

Yes 

Love 

c 

• 

- 

TABLE  IV. 3 
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RESOLUTE 


Trace  i 

ff  Per 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh . 

ss  -power 

49 

7.1 

.  81 

0.5 

ZN 

Yes 

No 

6.9 

0.79 

40 

EZ 

Yes 

No 

6.9 

0.70 

-35 

NE 

Yes 

No 

5.6 

0.56 

-55 

ZN 

Yes 

Yes 

4.3 

0.88 

-3 

NE 

Yes 

Yes 

Love 

50 

6.9 

0.74 

6.7 

NE 

Yes 

Yes 

Love 

6.6 

.34 

-104 

ZN 

Yes 

Yes 

Rayleigh 

4.3 

.  80 

2.0 

NE 

Yes 

Yes 

Love 

4.7 

189 

16 

EZ 

Yes 

Yes 

51 

6.4 

.52 

111 

EZ 

Yes 

Yes 

6.0 

.45 

25 

NE 

Yes 

Yes 

Love  (?) 

52 

4.1 

0.71 

-88 

EZ 

Yes 

No 

Rayleigh 

5.6 

.42 

-105 

EZ 

Yes 

Yes 

Rayleigh 

4.8 

0.71 

0.7 

NE 

Yes 

Yes 

Love 

6.9 

0.50 

3.4 

NE 

Yes 

Yes 

Love 

53 

6.2 

.47 

-83 

ZN 

Yes 

Yes 

Rayleigh 

6.4 

.41 

4.3 

NE 

Yes 

Yes 

Love 

4.7 

189 

-10 

NE 

Yes 

Yes 

54 

6.1 

.43 

-95 

ZN 

Yes 

Yes 

Rayleigh 

4.4 

189 

1.8 

NE 

Yes 

Yes 

Love 

6.4 

.53 

16 

NE 

Yes 

Yes 

Love 

55 

6.6 

.60 

96 

NE 

Yes 

Yes 

4.9 

187 

-104 

EZ 

Yes 

Yes 

Rayleigh 

• 

• 

race 

#  Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

ss  -power 

56 

4.7 

.  84 

-86 

ZN 

Yes 

Yes 

Love 

57 

6 . 6 

.63 

-5.3 

NE 

Yes 

Yes 

Love 

4.8 

.75 

93 

EZ 

Yes 

Yes 

6 . 6 

.  55 

46 

EZ 

Yes 

Yes 

58 

6.6 

.  80 

15 

EZ 

Yes 

Yes 

5.8 

.  80 

3.9 

NE 

Yes 

Yes 

Love 

59 

6.2 

.70 

-8 

EZ 

Yes 

Yes 

4.9 

.90 

0.5 

NE 

Yes 

Yes 

Love 

60 

4.3 

.92 

1.1 

NE 

Love 

« 

TABLE  IV. 4 
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OTTAWA 


Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh . 

ss  -power 

40 

7.3 

.60 

-110 

EZ 

Yes 

Yes 

Rayleigh 

41 

6.6 

.43 

-80 

EZ 

Yes 

Yes 

Rayleigh 

6.6 

.68 

49 

NE 

Yes 

Yes 

4.5 

.70 

-95 

EZ 

Yes 

Yes 

Rayleigh 

5.8 

.37 

2.6 

NE 

Yes 

Yes 

Love 

4.5 

.51 

2.5 

NE 

Yes 

Yes 

Love 

42 

6.0 

.64 

-37 

EZ 

Yes 

Yes 

6 . 6 

.  61 

52 

NE 

Yes 

Yes 

43 

5.7 

.42 

-72 

EZ 

Yes 

Yes 

Rayleigh 

7.1 

.  50 

-16 

NE 

Yes 

Yes 

-  ) 

Love 

4.2 

.59 

-1.2 

NE 

Yes 

Yes 

Love 

44  5 

6.6 

.56 

29 

EZ 

Yes 

Yes 

5.6 

.48 

75 

EZ 

Yes 

Yes 

Rayleigh 
(First  Mode) 

4.6 

.62 

-15 

NE 

Yes 

Yes 

Love 

45 

7.1 

.74 

96 

ZN 

Yes 

Yes 

7.1 

.57 

119 

EZ 

Yes 

Yes 

7.1 

.42 

134 

NE 

Yes 

Yes 

Love  (?) 

46 

6.6 

.51 

166 

EZ 

Yes 

Yes 

6.6 

165 

167 

ZN 

Yes 

Yes 

6 . 6 

.44 

8.4 

NE 

Yes 

Yes 

Love 

5.2 

.76 

-92 

EZ 

Yes 

Yes 

Rayleigh 

5.2 

.44 

156 

NE 

Yes 

Yes 

Love 

« 

4. 

- 
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Trace  # 

Per¬ 

Coher¬ 

Phase  Comp¬ 

Peak 

Peak  Cro¬ 

Type 

iod 

ency 

onent 

Coh. 

ss  -power 

47 

4.9 

.  66 

-100  EZ 

Yes 

Yes 

Rayleigh 

5.3 

.44 

-92  ZN 

Yes 

Yes 

Rayleigh  (?) 

6.9 

.60 

47  EZ 

Yes 

Yes 

6.9 

.57 

14  NE 

Yes 

Yes 

Love 

TABLE  IV. 5 
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EDMONTON 


Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  cro¬ 

Type 

iod 

ency 

onent 

Coh . 

ss  -power 

1 

7.1 

.  87 

-3.3 

NE 

No 

Yes 

Love 

3.8 

.92 

-3.0 

NE 

Yes 

Yes 

Love 

2 

3 

5.3 

.  87 

1.5 

NE 

Yes 

Yes 

Love 

7.1 

.  32 

-119 

EZ 

No 

Yes 

4 

7.1 

.53 

108 

EZ 

Yes 

Yes 

5 

4.0 

.74 

154 

NE 

Yes 

No 

4.0 

.96 

-173 

EZ 

Yes 

No 

6 

5.8 

.  55 

-100 

ZN 

Yes 

Yes 

6.4 

.45 

74 

NE 

Yes 

Yes 

Love  (?) 

7 

7.1 

0.49 

-83 

EZ 

Yes 

No 

Rayleigh 

8 

5.3 

.75 

1.9 

NE 

Yes 

Yes 

Love 

4.0 

.71 

-160 

EZ 

Yes 

Yes 

9 

3.2 

.98 

-3.1 

NE 

Yes 

Yes 

Love 

4.0 

.81 

164 

EZ 

Yes 

Yes 

10 

3.  8 

.87 

-174 

NE 

Yes 

Yes 

Love 

4.3 

.  89 

-11 

ZN 

Yes 

Yes 

11 

4.3 

.  69 

-112 

EZ 

No 

Yes 

Rayleigh 

5.8 

.60 

-43 

NE 

Yes 

Yes 

Love 

5.8 

.  35 

-109 

ZN 

Yes 

Yes 

Rayleigh 

♦ 

- 
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Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  cro¬ 

Type 

iod 

ency 

onent 

Coh. 

ss  -power 

12 

6.4 

.46 

-12 

NE 

Yes 

Yes 

Love 

6.4 

.  50 

33 

EZ 

Yes 

Yes 

13 

8.0 

.53 

-6.7 

NE 

Yes 

Yes 

Love 

14 

3.2 

.99 

0.9 

NE 

Yes 

Yes 

Love 

15 

16 

8.0 

.  82 

3 

NE 

No 

Yes 

Love 

17 

4.0 

.97 

-169 

EZ 

Yes 

Yes 

2.5 

.53 

-102 

EZ 

Yes 

Yes 

Rayleigh 

18 

5.3 

.59 

-15 

NE 

Yes 

Yes 

Love  (?) 

2.5 

.55 

-77 

EZ 

Yes 

Yes 

Rayleigh 

4.6 

.  89 

-164 

EZ 

Yes 

Yes 

19 

4.9 

.62 

-95 

EZ 

Yes 

Yes 

Rayleigh 

3.8 

.90 

8.6 

NE 

Yes 

Yes 

Love 

20 

6.4 

.43 

100 

EZ 

Yes 

¥es 

Rayleigh  (?) 
(First-Mode) 

21 

6.4 

.50 

-91 

EZ 

Yes 

Yes 

Rayleigh 

5.3 

.52 

163 

NE 

Yes 

Yes 

Love 

22 

7.1 

.62 

-117 

EZ 

Yes 

Yes 

Rayleigh 

7.1 

.77 

-13 

NE 

No 

Yes 

Love 

3.2 

.90 

1.9 

NE 

Yes 

Yes 

Love 

23 

6.4 

.49 

-13 

NE 

Yes 

Yes 

Love 

2.6 

.95 

-1.0 

NE 

Yes 

Yes 

Love 

K 
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Trace  # 

Per¬ 

Coher¬ 

Phase 

Comp¬ 

Peak 

Peak  cro- 

Type 

iod 

ency 

• 

onent 

Coh. 

fs  -power 

24 

5.8 

CO 

* 

-46 

EZ 

Yes 

Yes 

25 

6.4 

00 

LO 

• 

-70 

EZ 

Yes 

Yes 

Rayleigh 

28 

5.3 

.70 

17 

NE 

Yes 

Yes 

Love 

31 

5.8 

.36 

-63 

EZ 

Yes 

Yes 

Rayleigh 

5.8 

.38 

7 

NE 

Yes 

Yes 

Love 

32 

6.4 

.57 

13 

NE 

Yes 

Yes 

Love 

33 

6.4 

.81 

-3 

NE 

Yes 

Yes 

Love 

6.4 

.  36 

-85 

EZ 

Yes 

Yes 

Rayleigh 

38 

o 

• 

00 

.34 

-103 

EZ 

Yes 

Yes 

Rayleigh 

7.1 

.62 

3 

NE 

No 

Yes 

Love  (?) 

4 
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isfied.  An  event  with  a  peak  on  coherency  and  cross¬ 
power,  and  a  phase  of  0°  or  180°  +20°  between  two  hor¬ 
izontal  components  could  be  tentatively  identified  as  a 
Love  wave  if  no  vertical  peak  were  present.  First  mode  of 
Rayleigh  waves  could  be  identified  with  a  peak  on  coherency 
and  crosspower,  and  a  phase  of  +90°  +20°  between  the  vertical 
and  a  horizontal  component.  Since  all  the  three  criteria 
are  seldom  satisfied  simultaneously,  use  was  made  of  either 
two  or  one  of  the  criteria.  When  only  one  criterion  was 
used,  care  was  taken  to  see  that  the  other  two  were  nearly 
satisfied,  and  the  instances  where  only  one  criterion  was 
used  are  sparse. 

In  all  the  cases  where  a  Rayleigh  wave  was  identified 
the  coherency  between  the  vertical  and  a  horizontal  component 
rarely  exceeds  0.7  but  is  often  around  0.5  where  as  in  most 
cases  when  a  Love  was  identified  the  coherencies  generally 
exceed  0.7  and  more  often  they  are  around  0.9.  This  again 
shows  the  complex  nature  of  the  microseisms.  A  large  amount 
of  noise  on  the  vertical  component  could  be  the  reason  for 
such  low  coherency.  The  word  "noise"  has  to  be  understood 
here  as  any  energy  interfering  with  the  Rayleigh  wave  signal. 
A  body  wave  (for  example  a  long  period  leaking  P-mode)  reg¬ 
istered  on  the  vertical  component  with  the  same  period  would 
definitely  have  a  low  coherence  with  the  horizontal  component 
when  other  wave  types  are  present  at  the  same  time.  Theor¬ 
etically,  for  an  LP-mode  the  vertical  leads  the  horizontal 


. 
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component  by  180°  but  for  a  Rayleigh  mode  it  lags  behind 
the  horizontal  component  by  +90°.  So  depending  on  the 
path  difference,  the  phase  between  the  vertical  component 
and  the  horizontal  component  would  vary  between  these  two 
limits  if  the  signal  was  composed  of  LP-mode  and  Rayleigh 
mode  with  approximately  the  same  period.  Thus  the  low  coh¬ 
erency  between  these  two  components  or  a  lack  of  peak  on 
the  crosspower  between  them  or  a  wrong  phase  value  between 
them  is  neither  due  to  the  recording  nor  due  to  the  method 
of  calculation,  but  is  certainly  due  to  the  arriving  micro¬ 
seisms  themselves.  As  will  be  shown  by  the  studies  on  the 
synthetic  microseisms,  a  large  amount  of  white  noise  would 

be  required  to  lower  the  coherency  and  change  the  phase  to 

) 

the  observed  values.  In  addition  to  these  factors,  low 
values  of  coherency  between  the  vertical  and  horizontal  com¬ 
ponent  in  question  could  also  be  attributed  to  the  fact  that 
the  autopower  peaks  at  slightly  different  periods  on  these 
two  instruments.  The  latter  is  a  further  evidence  which  re¬ 
inforces  the  idea  that  body  waves  of  periods  around  6  sec 
are  indeed  present  in  the  microseisms.  Even  though  the  prob¬ 
lem  of  identification  of  a  wave  type  is  beset  with  tremendous 
amounts  of  difficulty,  and  because  of  this  not  all  wave  types 
could  be  properly  identified,  unmistakable  Rayleigh  and  Love 
wave  examples  are  not  entirely  lacking  in  the  data. 

IV.  2  Five  examples  of  Rayleigh  wave  can  be  cited  at 


. 
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Victoria  at  the  periods  of  5.8,  4.5,  6.0,  6.0,  and  4.9 
sec.  The  maximum  value  the  coherency  reaches  at  these 
periods  is  0.47.  Out  of  the  five,  three  are  the  events 
involving  east-west  and  vertical  components.  Two  cor- 
resond  to  the  recordings  on  north-south  component  and  the 
vertical.  This  suggests  that  the  Rayleigh  waves  are  arriv¬ 
ing  at  Victoria  from  at  least  two  directions  roughly  per¬ 
pendicular  to  each  other,  i.e.  from  north  and  west.  Excel¬ 
lent  examples  of  Love  waves  occur  at  6.1  and  4.5  secs  on 
the  record  #13,  at  6.1  and  4.6  secs  on  record  #14,  6.4  and 
4.6  secs  on  record  #15,  and  5.8  secs  on  record  #16.  The 
coherencies  are  high,  usually  greater  than  0.8  between  east 
and  north  with  phases  close  to  0°  or  180°.  There  are  cer¬ 
tain  cases  where  the  coherency  between  EW  and  vertical  is 
very  high  (.90)  but  the  phase  is  only  -0.3°.  This  wave  is 
tentatively  identified  as  some  kind  of  a  body  wave.  Few 
cases  of  first  mode  of  Rayleigh  wave  could  be  clearly  iden¬ 
tified.  Some  important  conclusions  can  be  drawn  from  this 
analysis . 

1.  Surface  waves  occur  at  two  periods,  one 
around  6.0  seconds  and  the  other  around  4.5 
seconds . 

2.  Rayleigh  and  Love  waves  arrive  at  the  station 
in  the  same  time  segment  having  the  same  fre¬ 
quency. 

3.  High  coherency  between  the  two  horizontals 
shows  the  Love  wave  is  recorded  almost  with¬ 
out  any  interference.  Low  coherency  between 
the  horizontal  and  vertical  indicates  some 
kind  of  interference  for  the  Rayleigh  type 


' 


vr 


61 


of  wave.  It  could  be  different  modes  of 
of  Raleigh  wave  or  in  addition  there  could 
be  ayf  iP-type  of  motion  as  the  phase  values 
suggest. 

4.  Direction  of  approach  is  not  very  clear  from 
this  analysis,  but  westerly  direction  seems 
to  be  predominant. 


IV. 3  Rayleigh  wave  with  periods  of  6.7,  6.2,  and  4.7 
seconds  with  coherencies  between  east-west  and  vertical  .48, 
.45,  and  .51  and  phases  -76,  -111,  and  -76  respectively  is 
recorded  at  Penticton.  At  this  station  Love  waves  have  per¬ 
iods  of  6.7,  7.1,  6.1,  4.6,  4.5,  6.4,  4.6,  6.4,  4.6,  6.0, 

4.5  seconds  with  high  coherencies  between  north-south  and 
east-west  components  and  right  order  of  magnitude  for  the 
phases.  The  coherencies  are  as  high  as  .93  and  phases  are 

Q 

very  approximately  either  zero  or  180  .  Thus  we  have: 

1.  Both  types  of  surface  waves  are  again  present 
at  two  periods,  one  ranging  between  6. 0-7.0 
and  the  other  around  4.5  seconds. 

2.  They  occur  in  the  same  time  segment. 

3.  There  is  a  low  coherency  between  vertical  and 
horizontal  for  Rayleigh  waves  and  high  coher¬ 
ency  between  the  two  horizontals  for  Love  waves. 
Phase  values  for  Rayleigh  are  out  by  -15  to 
-21°;  phase  values  for  Love  are  of  right  order. 

4.  Direction  of  approach  seems  to  be  somewhere 
in  the  north-west  quadrant. 

IV. 4  The  best  example  of  Rayleigh  wave  at  Resolute  occurs 
at  4.7  secs  with  a  coherency  of  .84  between  vertical  and 
north-south  with  a  phase  difference  of  -86  .  Another  example 


is  that  of  a  4.1  sec  event  with  a  coherency  of  0.71  and 
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phase  of  -88°  for  the  east-west  and  vertical.  Rayleigh 
wave  with  a  period  of  6.2  sec  with  a  coherency  of  .47  and 
phase  of  -83°  between  north-south  and  vertical  is  also 
recorded.  There  is  however  no  evidence  for  a  6.0  second 
Rayleigh  wave  coming  from  west  or  east.  Love  waves  of  6.9, 
6.4,  6.1,  and  4.5  seconds  are  also  recorded  at  Resolute. 

They  have  expected  phase  relationships  between  the  two 
horizontals . 

IV. 5  Two  events,  one  around  6.5  sec  and  the  other  around 
4.6  sec  are  found  on  Ottawa  records.  Rayleigh  waves  have 
an  east-west  direction  of  approach.  One  excellent  example 
occurs  on  record  #46  at  the  period  of  5.2  sec  with  a  coher- 
ency  of  0.76  and  phase  of  -92°  between  the  vertical  and  east- 
west.  Love  and  Rayleigh  waves  have  approximately  the  same 
order  of  coherency.  There  is  one  instance  of  a  northerly 
approaching  Rayleigh  wave. 

IV.  6  The  best  example  of  Rayleigh  wave  occurrence  at 
Edmonton  is  on  record  #21  at  period  of  6.4  sec  with  a  coh¬ 
erency  of  .50  between  east-west  and  vertical  components  and 
a  phase  of  -91°.  In  addition  to  the  event  with  a  period  of 
4.9  sec.,  another  event  at  2.5  sec  can  be  easily  seen  at  this 
station.  Surface  wave  period  changes  with  time  and  the  range 
of  variation  is  5. 3-8.0  sec  and  4. 0-5.0  sec  and  the  third 
one  occurring  around  2.5  sec.  This  wide  range  is  probably 
due  to  the  crustal  layering  at  Edmonton.  The  2.5  sec  event 
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is  likely  due  to  the  ^consolidated  sediments  overlying 
the  precambrian.  The  relatively  larger  amplitudes  of  the 
short  period  peak  at  Edmonton  compared  to  stations  at  both 


greater  and  lesser  distances  from  the  presumed  source  sug¬ 
gests  that  the  effect  cannot  be  due  solely  to  differing 
attenuations  with  distance  but  must  be  related  in  some  way 
to  the  ^consolidated  sediments  overlying  the  Precambrian 
at  Edmonton. 

In  summary  then: 

1.  Rayleigh  and  Love  waves  arrive  at  the  station 
simultaneously. 

2.  They  have  approximately  the  same  periods. 

3.  Two  predominant  periods  can  be  recognized  for 
these  waves.  One  train  of  waves  has  a  period 
around  6.5  sec  and  the  other  around  4.5  sec. 

Yet  another  event  of  2.5  sec  could  be  identified 
only  at  Edmonton. 

4.  There  seems  to  be  more  than  one  direction  from 
which  the  waves  are  approaching  the  recording 
site. 

5.  The  coherency  between  the  two  horizontals  is 
high  indicating  a  coherent  propagation  of 

Love  waves  in  the  crust.  The  coherency  between 
the  vertical  and  any  one  horizontal  component 
rarely  exceeds  0.70  suggesting  an  interference 
of  fundamental  and  higher  Rayleigh  modes  and 
AP-modes . 

6.  The  phases  are  in  most  cases  correct  for  a 
Love  wave  but  vary  by  approximately  20  for 
the  Rayleigh  wave,  a  circumstance  supporting 
the  interference  above. 

Certain  wave  types  with  a  coherency  of  as  high 
as  .90  between  the  vertical  and  a  horizontal 
component  but  phase  of  0°  or  180°  were  cons¬ 
picuous  in  the  spectral  analysis.  These  wave 
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types  have  a  period  also  around  6.5  sec.  and 
strongly  suggest  the  existence  of  long  per¬ 
iod  body  waves  in  microseisms. 

8.  One  or  two  examples  of  first  mode  of  Rayleigh 
wave  could  be  identified. 

IV. 7  Wherever  there  was  a  good  example  of  Rayleigh  wave, 
the  vertical  to  radial  ratio  was  calculated  (radial  being 
the  EW  component)  for  both  events  of  6.5  secs  and  4.5  secs 
and  the  results  are  listed  in  Table  IV. 6.  The  crustal  struc¬ 
ture  at  all  these  stations  is  approximately  known  from  the 
literature.  From  the  known  crustal  structure  Rayleigh  wave 
dispersion  curves  were  calculated  and  the  vertical  to  horiz¬ 
ontal  ratio  computed.  These  theoretically  calculated  ratios 
are  compared  with  the  observed  ones  in  Table  IV. 7.  The 
crustal  models  used  are  shown  in  Table  IV. 8.  The  models  at 
Victoria  and  Penticton  are  taken  from  White  et  al  (1965) , 
that  at  Edmonton  from  Cumming  et  al  (1966) ,  at  Ottawa  from 
Brune  et  al  (1963)  and  finally  at  Resolute  from  Sander  et 
al  (1965) .  The  method  used  in  computing  the  dispersion 
curves  is  that  described  by  Thomson  (1950)  and  later  gener¬ 
alized  by  Haskell  (1953)  and  subsequently  known  as  the  Haskell- 
Thomson  method.  A  program  earlier  described  by  Press  et  al 
(1961)  and  later  expanded  by  Harkrider  and  Anderson  (1962) 
was  used  to  calculate  the  dispersion.  An  isotropic  model 

was  used  for  all  locations. 

For  the  longer  period  event  the  agreement  between 


■ 
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VERTICAL-HORIZONTAL  AMPLITUDE  RATIOS  FOR  RAYLEIGH  WAVE 


Site 

Vertical2 

2 

Horizontal 

v2/h 

2 

V/H 

Period  in  Secs 

6.5 

4.5 

6.5 

4.5 

6 . 5 

4.5 

6.5 

4.5 

Victoria 

8.6 

.55 

15.2 

2.2 

.57 

.25 

.  75 

.5 

Penticton 

6 . 6 

5.5 

2.5 

6.0 

2.6 

.92 

1.6 

.96 

Edmonton 

3.1 

.74 

2.7 

1.8 

1.1 

.41 

.10 

.64 

Resolute 

.09 

.096 

.10 

.11 

.9 

.  87 

.95 

.93 

Ottawa 

.05 

.20 

.09 

.08 

.56 

215 

.  75 

1.6 

TABLE  IV. 7 


COMPARISON  WITH  THEORETICAL  VALUES 


Site 

6.5 

sec 

4.5 

sec 

Exper 

Theo 

Exper 

Theo 

Victoria 

0.75 

0.90 

0.50 

0. 80 

Penticton 

1.6 

0.98 

0.96 

0.98 

Edmonton 

1.01 

1.1 

0.64 

0.91 

Resolute 

0.95 

0.69 

0.93 

0.70 

Ottawa 

0.75 

0.72 

1.6 

0.66 

• 

' 

TABLE  IV. 8 
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CRUSTAL  MODELS 


1.  Victoria 


2.  Penticton 


Thickness 
2 . 5  km 
6 . 2  km 
45.0  km 


P-wave  Velocity 
4 . 4  km/sec 
6.0  km/sec 
6.73  km/sec 


7.74  km/sec 
1  km  sediments 


20  km 

5.91 

.  km/sec 

5  km 

6.8 

km/sec 

7.  8 

km/sec 

3. 

Edmonton 

2.0 

km 

3.0 

km/sec 

3.0 

km 

5.2 

km/sec 

10.0 

km 

6.1 

km/sec 

20.0 

km 

6.5 

km/sec 

11.0 

km 

7.2 

km/sec 

8.2 

km/sec 

4. 

Resolute 

20.0 

km 

6.0 

km/sec 

20.0 

km 

7.3 

km/sec 

8.2 

km/sec 

5. 

Ottawa 

6.0 

km 

5.64 

km/sec 

10.5 

km 

6.15 

km/sec 

18.7 

km 

6.6 

km/sec 

8.1  km/sec 


Comments 


1.  Victoria:  Travel  time  curves  of  White  and  Savage  (op.cit.) 
all  indicate  that  a  thin  layer  of  low  velocity  material 
is  present.  Values  used  are  estimated  from  this  data. 


2.  Penticton:  Negligible  low  velocity  surface  material  is 
present.  (White,  Bone  and  Milne,  1968). 

3.  Edmonton:  This  model  is  based  on  Cumrning  and  Kanasewich 
(op.cit.)  Ellis  and  Basham  (op.  cit.)  indicate  thinner 
second  layer. 
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observed  and  calculated  values  is  good  at  all  stations 
except  at  Penticton.  For  the  shorter  period  event  the 
agreement  is  good  only  at  Penticton  and  Resolute.  Major 
discrepancies  occur  for  the  long  period  event  at  Penticton. 

For  the  shorter  period  event  the  agreement  is  good  only 
at  Penticton  and  Resolute.  Major  discrepancies  occur  for 
the  long  period  event  at  Penticton  and  short  period  event 
at  Ottawa.  If  we  assume  microseisms  contain  body  waves  and 
surface  waves,  this  discrepancy  could  be  attributed  to  a.;, 
major  contribution  from  ^P-waves  which  would  reflect  in  the 
larger  value  of  amplitude  on  the  vertical  component.  Observe 
that  the  values  cited  in  Table  IV. 6  pertain  to  only  one 
record  where  Rayleigh  wave  was  identified  on  the  basis  of 
correct  phase  and  coherency  values.  Discrepancy  could  then 
be  ascribed  to  the  lack  of  knowledge  of  crustal  layering  at 
Penticton.  However  the  discrepancy  at  Ottawa  for  the  4.5 
sec  event  can  not  be  readily  explained.  These  results  at 
least  partly  confirm  that  the  radial  direction  is  approxr 
imately  EW  at  every  station  for  the  record  chosen  with  the 
waves  approaching  the  sites  from  a  westerly  direction. 

IV. 8  All  the  results  enumerated  so  far  point  out  that  the 
arriving  microseisms  are  very  complex  in  their  nature.  They 
arrive  from  different  directions;  they  are  not  simple  Rayleigh 
and  Love  waves,  but  are  a  mixture  of  fundamental  and  higner 
modes;  the  wave  types  even  include  some  body  waves.  Since 
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the  coherency  is  usually  less  than  0.7  (especially  for 
Rayleigh  waves) ,  there  appears  to  be  considerable  random 
noise  contained  in  the  records.  It  is  this  randomness 
that  seems  to  be  the  crux  of  the  problem. 

An  attempt  was  made  to  combine  a  signal  with  random 
numbers  and  calculate  the  crosspower,  coherency  and  phase 
between  two  components  simulating  a  Rayleigh  wave,  with  -90° 
phase  difference,  and  between  a  second  set  of  two  components 
with  180°  phase  difference  simulating  a  Love-type  of  wave. 
Random  numbers  were  generated  and  added  to  these  signals  in 
different  proportions.  Five  cases  were  tested. 

The  signal  to  (white)  noise  ratio  in  these  five  cases 
is  1)  1/1,  2)  1 /J  2,  3)  1/5,  4)  1/10,  and  5)  1/20. 

Table  IV. 9  lists  the  phase  and  coherency  values  between  the 
two  traces  for  these  five  cases  for  the  two  wave  types  con¬ 
sidered.  The  following  features  stand  out  of  these  investig¬ 
ations  : 


1.  Coherency  falls  to  less  than  one  half  for  both 
wave  types  only  in  the  case  when  the  noise  is 
20  times  the  signal. 

2.  Phase  changes  by  10°  for  the  Rayleigh  type 

of  motion  and  16°  for  the  Love  type  of  motion 
in  the  fifth  case  when  the  noise  is  20  times 
the  signal. 

3.  The  values  obtained  for  coherency  in  the  fifth 
case  (a  signal  to  white  noise  ratio  of  1/20) 
are  comparable  to  the  experimentally  obtained 
values  for  the  Rayleigh  type  of  motion. 

The  comparison  of  the  experimental  values  with 
the  ones  obtained  in  Table  IV. 9  for  coherency 
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TABLE  IV. 9 
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SYNTHETIC  MICROSEISMS  AND 

RANDOM  NOISE 

S/N 

Rayleigh 

Type 

Love  Type 

Coherency 

Phase 

Coherency 

Phase 

1 

1.0 

-89 

1.0 

179.0 

1  /si  2 

0.99 

i 

00 

00 

• 

o 

.99 

178.4 

1/5 

0.84 

o 

• 

00 

l 

.90 

176.0 

1/10 

0.78 

l 

00 

to 

• 

o 

.  86 

170.0 

1/20 

0.44 

-80.9 

.49 

164.8 
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of  Love  type  motion  indicates  that  in  the 
case  of  Love  wave  the  signal  to  noise  ratio 
is  somewhere  between  1/n /  2  and  1/5. 

5.  The  comparison  of  phase  values  would  put  the 
signal  to  noise  ratio  even  lower  for  both 
wave  types.  Even  a  considerable  amount  of 
white  noise  by  itself  is  not  capable  of  chang¬ 
ing  the  phase  by  more  than  20°. 

6.  The  microseisms  as  observed  on  the  records 
exhibit  a  clear  sinusoidal  motion  with  beat 
phenomenon  superimposed  on  them.  A  signal 
to  noise  ratio  of  1/20  can  not  produce  the 
type  of  sinusoids  seen  on  the  records  and 
hence  this  type  of  noise  clearly  does  not 
model  the  observed  records. 

7.  In  order  to  explain  the  low  coherency  observed 
for  Rayleigh  waves,  and  at  the  same  time  the 
high  coherency  observed  for  Love  waves ,  one 
has  to  concede  that  white  noise  in  any  prop¬ 
ortion  is  not  an  adequate  explanation  of  the 
observed  data.  The  phase  values  suggest  that 
other  types  of  motion  must  be  present  in  the 
microseisms  but  associated  with  Rayleigh  and 
Love  waves.  That  these  spurious  phase  values 
are  not  introduced  by  the  method  of  analysis 
is  clear  from  the  fact  that  even  though  the 
data  on  all  the  components  are  subjected  to 
the  same  mathematical  operations,  the  phase 
between  horizontal  components  is  often  approp- 
iate  for  Love  waves,  while  results  including 
vertical  motion  are  seldom  reasonable  for 
simple  Rayleigh  waves. 

8.  Since  the  "randomness"  or  "noise"  is  not  solely 
due  to  1)  white  noise  which  is  incapable  of 
producing  the  observed  coherency  and  phase 
values  while  preserving  the  general  appearance 
of  the  records,  2)  the  method  of  analysis,  and 
3)  the  recording  instruments;  it  must  be  assoc¬ 
iated  with  the  nature  of  the  arriving  micro¬ 
seisms  themselves. 

9.  In  tables  IV. 1-5,  one  can  observe  some  wave 
types  which  have  a  phase  relationship  of  0 
or  180°  between  the  vertical  and  horizontal 
components  with  a  coherency  as  high  as  0.9. 

If  this  is  some  kind  of  a  leaking  LP-mode, 
then  the  interference  of  Rayleigh  motion  and 
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LP-motion  may  be  responsible  for  the  low 
coherencies  observed  between  the  vertical 
and  horizontal  components  for  the  majority 
of  the  data. 


. 


CHAPTER  V 


JENSEN'S  DIAGRAMS  AND  PARTICLE  MOTION  TRAJECTORIES 

V. 1  A  brief  description  of  the  method  of  plotting  Jensen's 
diagrams  was  given  in  1.6.  The  first  derivatives  of  the 
two  horizontal  components  with  respect  to  time  are  cal¬ 
culated  at  the  points  of  maxima  on  the  vertical  component. 
These  two  are  then  plotted  against  each  other  and  this 
plot  would  then  describe  the  velocity  of  the  particle. 

An  average  of  several  velocity  vectors  would  give  the 
direction  of  approach.  All  the  points  would  lie  on  the 
half  plane  away  from  the  source  for  a  fundamental  Rayleigh 

i 

mode,  because  of  the  retrograde  nature  of  the  particle 
motion.  For  a  first  Rayleigh  mode  they  would  lie  on  a 
half  plane  towards  the  source.  All  these  points  tend  to 
cluster  around  the  line  along  which  the  pure  Rayleigh  wave 
is  propagated.  For  a  combination  of  fundamental  Rayleigh 
mode  and  Love  wave  the  points  still  fall  on  a  half  plane 
away  from  the  source  but  in  this  case  the  points  are  widely 
dispersed  from  the  line  representing  the  direction  of 
approach.  For  a  combination  of  first  Rayleigh  mode  and 
Love  wave  the  points  still  lie  on  a  half  plane  towards  the 
source  but  they  are  widely  dispersed  around  the  line  along 
which  the  waves  are  propagated.  For  an  SV-type  or  SH-type 
of  motion  these  points  representing  the  velocity  vectors 
would  fall  on  the  half  plane  towards  the  source,  but 


' 
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form  a  straight  line.  In  the  former  case  there  is  no 
appreciable  horizontal  motion  and  in  the  latter  case  there 
is  no  appreciable  vertical  motion. 

If  the  points  are  spread  around  the  centre  (point 
of  intersection  of  the  two  axes  representing  slopes  of 
two  horizontal  components)  it  would  indicate  i)  either  the 
waves  are  approaching  the  site  from  different  directions, 
or  ii)  the  signal  is  composed  of  different  wave  types. 

It  is  desirable  to  suppress  the  noise  as  much  as 
possible  in  order  to  draw  any  meaningful  Jensen's  diagrams 
so  a  Butterworth  band  pass  filter  was  used  on  the  data 
before  plotting.  The  band  pass  was  centered  at  around  6.0 
sec  which  is  the  most  prevalent  period  in  the  microseisms. 
The  filter  designed  by  Butterworth  (1930)  is  fully  des¬ 
cribed  by  Weinberg  (1962) .  A  brief  outline  is  given  here. 

V. 2  The  function 

1 

I  Y  (^)  I  2=  1  +  Q  2n 

describes  the  low  pass  Butterworth  filter,  where  ^  is 
normalized  frequency,  ^  =  to/  being  the  cut  off  of 

frequency.  The  function  does  not  tend  to  oscillate  in  the 
band  pass  and  this  property  is  exploited  to  design  a  max¬ 
imally  flat  curve,  which  falls  off  rapidly  outside  the  pas 
band  as  n  is  increased.  Poles  of  the  low  pass  filter  are 
obtained  by  taking  the  [Laplace  transform  of  (1)  . 


' 


74 


Thus  | Y  (p)  | 


2  _ 


1  +  (-1) n  p2n 


(2) 


where  p  — p  +  iQ  and  p  =  o.  There  are  no  zeroes  in  a  low 
pass  filter  except  at  infinity  and  the  poles  are  given  by 

1  +  (-1)  np2n  =  0 


In  general  |Y^(p)|2can  be  represented  as 


(3) 


[YL(p)] 


Bn(p)B*  (p) 


B  (p)  being  known  as  Butterworth  polynomial, 
p  plane  B  (p)  can  be  shown  to  be 


In  the 


(4) 


Bn(P} 


n-1 

n  (p-e 
v=0 


ie(2v  +  1  +  n)) 


(5) 


For  a  high  pass  filter  the  frequency  transformation  is 
given  by  S  =3  +  iw  where  w  =  u)/u)t-,a)t,  is  cut  off  frequency. 

n  n 

For  a  band  pass  filter  the  normalized  frequency  transform¬ 
ation  is  p  =  S  +  i  (6) 

u 

The  normalized  relations  must  be  generalized  to  the  desired 
band  pass  filter  with  unnormalized  upper  and  lower  fre¬ 
quency  limits  co^  and  The  geometric  mean  is  sometimes 

called  the  center  frequency.  Any  other  frequency  with  a 
particular  attenuation  is  related  to  a  higher  frequency. 


by  the  relation 


Wo  =  /  0)k  " 

The  band  width  Aw  is  obviously  .  The  band  width 

will  be  rescaled  by  dividing  by  the  geometric  mean 


(7) 


6.  £  .  bo  •  •  •.  T  ;  . 


0)1 
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ft  i 


0)2  - 


0) 


Any  frequency  in  general  can  be  given  by 

ft  =  “  - 

0)  0  0) 

Normalized  frequency  is  then 

^ 2  ~  Q) 1 0)2 
ft  1  0)  (  0)2  “  0)1  ) 

Since  p  =  ift  and  S  =  iw,  we  obtain 

2 

S  -  P  ( 0)  2  ~  0)  1  )  S  +  0)  1 0)  2  =  0  (8) 

is  a  quadratic  in  S. 

The  transfer  function  has  the  form 

Vp)  “  (p-p1)  (p-p2)  (p-p3)  (p-p4)  (9) 

for  a  fourth  order  Butterworth  polynomial,  if  we  consider 
a  band  pass  filter,  p  =  S  +  Equation  (9)  can  be 

reduced  to 

YBP^S)  =  (S2b1S+c1^  (Sz+b2S+c2)  (Sz+b3S+c3)  (S2+b4S+c4) 

The  coefficients  b.  and  c.  are  obtained  from  the  pole 

J  J 

positions,  which  occur  in  complex  conjugate  pairs  to  make 
each  of  the  above  terms  in  brackets  real. 

(S  -  S.)  (S  -  S.*)  =  S2  +  b.S  +  c.  (11) 

J  D  D  D 

The  continuous  S  plane  representation  has  now  to  be  modi¬ 
fied  to  a  digital  Z  -  plane.  The  standard  Z  -  transform  is 


(10) 
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is 


Z  =  e-S(At) 


(12) 


(At)  being  the  digitizing  interval.  In  order  to  eliminate 
the  aliasing  errors  present  in  the  standard  z  -  transform, 
one  uses  a  bilinear  Z  -  transform  following  Golden  and 
Kaiser  (1964)  .  The  desired  transformation  is  then 

Sd  =  (2/ At)  (1-Z ) /  (1+Z )  (13) 

Sd  =  2/At  tanh  l/2At  because  of  (12) 

In  the  frequency  domain  Sd  =  iood  and  using  tanhx  =  -i  tanix. 


we  have  co.  =  jt- 
d  At 

From  (10)  then 
S  . 


tan 


At 


(14) 


~a  -  z2) 


Sj+bTsT+cT 

d  j  d  j 


1  +  Z  ,c  .  At  ~  4  x  Z  2  , 2 

—  (  3  — )  +  —  ( 


At' 


b  .  .  c  .  Atx 
a.  'At  -  3  + 


where 


= 

i.d 

3 

- 

Z2) 

2 

+  b. 
3 

B  . 

3 

At 

aj  ~  At 

+ 

c .— =■ 
3  2 

B.  =  1- 
3 

D-,  .Z 

13 

+ 

D0  .  Z 
23 

(15) 


Z  -  transform  of  a  Butterworth  band  pass  filter  with  eight 
poles  is  given  by 

(1-Z2)4 

W(Z)  “  (Z)B2  (Z)B3  (Z)B4  (Z) 


(16) 


The  Z  -  transform  thus  has  both  singularities  and  poles. 
Equation  (16)  can  be  modified  to  give  a  Z  -  transform  for  a 

band  pass  filter  with  four  poles. 

2  2 

(1-zV 


(Z)B2  (Z) 


W  (Z) 


(17) 
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where  B1  (Z)  =  1  -  1.307476Z  +  0.450190Z2 
B2  (Z)  =  1  -  1.494986Z  +  0.724783Z2 
In  equation  (10) 

h±  =  94.10519;  c  =  2882.33984 
fc>2  =  40.34528;  c2  =  3975.04687 
The  singularities  in  equation  (17)  are  given  by 
Z1  =  1.452137  ±  i  0.335533 
Z2  =  1.031333  ±  i  0.562206 

A  programme  was  in  operation  previously  (Alpaslan, 
M.Sc  Thesis  1968)  for  an  eight  pole  transfer  function  for 
a  ButterWarth  band  pass  filter.  The  filter  was  thought  to 
be  too  severe  for  microseisms  and  the  programme  was  modi¬ 
fied  for  a  four  pole  transfer  function.  The  programme 
can  be  manipulated  to  give  a  zero  phase-shift  response. 

The  amplitude  response  is  fairly  flat  with  sharp  cut  offs 
at  the  two  limits  of  frequency. 

The  operation  of  this  filter  on  microseisms  is 
illustrated  in  Figures  V.l-3.  In  Figure  V.l  is  shown  the 
autopower  spectrum  of  vertical  component  of  a  particular 
record  (from  Ottawa)  on  which  two  peaks  at  around  6  sec¬ 
onds  and  16  seconds  were  found.  The  autopower  specturm 
was  calculated  for  this  trace  after  it  was  band  pass  fil¬ 
tered  for  two  bands  .10  -  .05  and  .20  -  .10,  and  the  re¬ 
sults  are  shown  in  Figures  V.2  and  V.3.  The  filter  is 
obviously  very  effective  in  rejecting  signal  outside  the 
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POWER  SPECTRUM  OF  THE  FILTERED  DATA 
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FIG  V. 3 
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desired  pass  band. 

V.3  The  microseismic  data  were  Butterworth  band  pass  fil¬ 
tered  and  Jensen's  diagrams  were  drawn  for  several  cases. 
These  diagrams  did  not  show  any  encouraging  results.  Since 
the  filter  did  not  introduce  any  undesirable  frequencies 
or  phases  (note  that  the  filter  was  manipulated  to  give  a 
zero  phase-shift  response) ,  the  lack  of  any  coherent  direc¬ 
tion  of  approach  is  attributed  to  the  arriving  microseisms 
themselves.  In  Figures  V.4-5  are  shown  some  typical 
Jensen's  diagrams  for  the  three  stations  Victoria,  Pen¬ 
ticton  and  Edmonton.  In  almost  all  these  diagrams  the 
points  are  not  confined  to  one  half  plane,  thus  indicating 
that  the  microseisms  recorded  at  these  sites  are  indeed 
composed  of  different  types  of  surface  waves  superimposed 
in  a  complex  fashion.  As  will  be  shown  in  the  next  chapter 
on  the  rotated  crosspower,  it  is  possible  to  find  in 
several  cases  a  direction  of  approach  correct  to  plus  or 
minus  ten  degrees.  Since  this  magnitude  of  error  is  not 
solely  capable  of  producing  the  large  amounts  of  scattering 
observed  in  the  Jensen's  diagrams,  it  is  believed  that 
microseisms  are  very  complex  in  their  content  of  physical 
wave  types.  If  it  were  possible  to  separate  the  different 
wave  types,  perhaps  using  an  array  of  tripartite  stations, 
then  meaningful  diagrams  might  be  drawn,  but  with  the 
present  data  it  appears  that  this  technique  is  not  capable 
of  producing  useful  results. 
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JENSEN’S  DIAGRAMS  AT  VICTORIA 
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JENSEN'S  DIAGRAMS  AT  PENTICTON 
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V. 4  It  is  also  possible  to  draw  particle  motion  trajec¬ 
tories  using  the  filter  traces.  Two  such  sets  of  diagrams 
can  be  drawn;  one  of  them  between  the  vertical  component 
and  the  horizontal  component  rotated  into  the  approximate 
direction  of  approach  of  the  waves,  and  the  second  between 
the  two  horizontal  components  again  rotated  so  that  the 
radial  direction  corresponds  to  the  approximate  direction 
of  source.  The  former  should  be  a  Rayleigh  wave  trajectory 
and  the  latter  a  Love  wave  trajectory  with  respectively 
elliptical  and  linear  particle  motions. 

In  figures  V.6-7  are  illustrated  the  particle  tra¬ 
jectories  drawn  between 

1)  the  vertical  and  the  radial  components  (for 
Rayleigh  wave) 

2)  the  radial  and  the  transverse  components  (for 
Love  wave)  at  Victoria  and  Penticton.  The  examples  shown 
were  chosen  from  among  several  such  diagrams.  Although 
particularly  in  Fig.  V. 7  there  is  some  indication  of  polari¬ 
zation  of  the  trajectories,  if  one  were  to  examine  similar 
diagrams  constructed  at  slightly  different  times  the  plane 
of  polarization  would  be  seen  to  have  rotated  to  a  new 
direction  and  in  general  there  does  not  seem  to  be  any 

clear  cut  particle  motion. 

It  may  also  be  noted  that  for  the  vertical-radial 

plot 

1)  the  axis  of  the  trajectory  is  inclined  to  the 
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PARTICLE  TRAJECTORIES  AT  VICTORIA 
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PARTICLE  TRAJECTORIES  AT  PENTICTON 


FIG  V.  8 


80 


vertical ; 

2)  the  motion  is  retrograde  at  the  start,  but  it 
switches  to  prograde  later  in  this  particular 
time  segment;  this  can  be  especially  seen  in 
Fig.  V.8. 

3)  there  is  a  large  ''vertical"  to  "horizontal" 
ratio  (almost  three  to  one) 


and 


% 


CHAPTER  VI 


DIRECTIONAL  SPECTRA 


VI. 1  It  is  clear  from  what  has  been  said  so  far  that 

the  determination  of  such  a  simple  property  as  the  direction 
of  approach  of  the  surface  waves  is  a  difficult  problem 
when  one  has  available  only  single  three-component  stations. 
It  is  even  more  difficult  in  most  cases  to  determine  with 
any  degree  of  certainty  what  types  of  waves  are  present 
in  the  data.  In  some  cases  one  can  assume  that  the  source 
is  roughly  in  the  direction  of  one  of  the  horizontal  com¬ 
ponents  and  thus  get  some  idea  of  the  possible  wave  types 
from  the  relations  between.1'  crosspower,  coherency  and  phase 
but  a  more  detailed  analysis  is  obviously  desirable. 

The  approach  taken  here  is  to  first  compute  the 
resultant  horizontal  motions  at  various  azimuths  starting 
in  a  westerly  direction  and  rotating  the  axes  clockwise 
at  10°  intervals  to  a  maximum  of  180°  at  which  point  the 
data  then  repeats  with  the  phase  reversed.  For  convenience 
the  initial  west  component  is  referred  to  as  "radial"  and 
the  orthogonal  component  as  "transverse",  and  this  term¬ 
inology  is  maintained  throughout  the  rotation. 

The  data  were  then  used  to  calculate  the  crosspower, 
coherency  and  phase  between  the  vertical  and  the  radial  com¬ 
ponents,  and  between  the  radial  and  the  transverse  com¬ 
ponents,  for  each  rotation.  These  quantities  were  written 
on  a  grid  as  functions  of  two  variables,  azimuth  and  period 
and  then  contoured.  Four  such  maps,  two  for  crosspower  and 
two  for  coherency  were  prepared  at  each  station.  In  addition 
to  coherency  and  crosspower,  the  phase  between  the  two  com¬ 
ponents  involved  was  also  written  on  the  grid  as  a  function 
of  two  variables,  azimuth  and  period.  No  contour  diagrams 
for  phase  values  was  drawn  as  they  varied  in  an  extremely 
random  fashion  from  period  to  period. 
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VI  •  2  Some  properties  of  these  contour  maps  can  be 

enumerated  as  follows: 

1.  Maxima  and  minima  in  the  crosspower 
and  coherency  of  the  original  data  will 
now  be  reflected  in  the  contour  map  as 
high  and  low  closures. 

2.  From  the  properties  of  the  rotation  and 
the  crosspower  calculations,  it  can  be 
seen  that,  when  dealing  with  vertical- 
horizontal  power,  the  value  will  repeat 
after  180  (one  full  rotation)  and  any 
high  will  have  a  corresponding  low  sep¬ 
arated  from  it  by  90°  (one  half  rotation) . 

3.  The  values  of  crosspower  and  coherency 
between  the  two  resolved  horizontal  com¬ 
ponents  will  be  repeated  every  90°  and 
the  highs  and  lows  will  be  separated  by 
45  . 

4.  For  vertical-radial  data  a  single  direction 
of  propagation  of  the  wave  will  result 

in  a  high,  sharply  defined  in  azimuth, 
whereas  a  diffuse  source  will  result  in 
a  broad  peak  and  the  half-power  width 
can  be  used  as  a  measure  of  the  angular 
size  of  the  source  region. 

5.  For  radial-transverse  power  a  Rayleigh 
wave  coming  from  a  single  direction  will 
result  in  a  sharply  defined  high  as  abl?vt, 
but  the  maximum  will  occur  45°  away  from 
the  direction  of  approach.  The  presence 
of  Love  waves  in  the  data  will  result  in 
departures  from  this  relationship  and  the 
position  of  maxima  and  minima  will  be  a 

a  function  of  the  magnitude  and  direction 
of  the  Love  component.  If,  for  example, 
the  Love  and  Rayleigh  waves  approach  from 
the  same  direction  and  are  coherent  the 
maximum  in  the  horizontal  crosspower  will 
occur  at  the  same  direction  as  that  for 
vertical-radial  crosspower  and  the  coh¬ 
erency  will  be  high. 

6.  The  crosspower  90°  from  a  maximum  represents 
noise  from  that  azimuth  and  if  we  make 

the  reasonable  assumption  that  the  noise 
is  isotropic  we  can  use  the  ratio  between 
high  and  low  crosspower  as  a  measure  of 
microseismic  signal  to  noise  ratio. 
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7.  The  ratio  of  the  magnitude  (A  )  of  the 

vertical  autopower  spectrum  to  the  cross¬ 
power  (AB)  between  vertical  and  radial  at 
the  same  period  is  the  ratio  of  vertical 
to  radial  amplitudes  and  this  may  then 
be  compared  with  Rayleigh  wave  values  cal- 
ulated  from  dispersion  curves  for  layered 
media. 

VI.  3  A  sample  of  these  contour  maps  for  each  station  is 

shown  in  Figures  VI. 1-5.  It  should  be  mentioned  here  that 
the  samples  chosen  correspond  to  those  times  at  which  there 
is  a  maximum  peak  value  on  the  autopower  spectrum  of  the 
vertical.  This  maximum  is  reached  at  different  times  at 
different  sites.  As  a  result  these  samples  also  differ  in 
time  and  the  position  of  the  storm  center  is  therefore  dif¬ 
ferent  for  the  various  examples.  On  these  grids  the  angles 
of  rotation  are  given  in  columns  and  the  periods  as  rows. 

By  the  side  of  the  contour  map  is  given  a  plot  of  the  auto¬ 
power  spectrum  of  the  vertical  component.  This  plot  is 
given  only  on  those  maps  where  vertical  component  is  involved 
either  in  the  crosspower  or  in  the  coherency  calculation. 

At  Victoria  there  is  one  peak  at  6.56  seconds  on 
the  autopower  spectrum  of  the  vertical  component.  A  high 
at  this  period  can  be  seen  on  both  contour  maps,  for  coh¬ 
erency  and  crosspower  between  the  vertical  and  the  radial 
components  at  the  directions  of  60  and  65  measured  clock¬ 
wise  from  West.  (Note  that  all  directions  mentioned  here 
are  measured  in  this  manner.)  Phase  for  this  event  is  90 
a  possible  indication  of  prograde  Rayleigh  motion.  On 
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these  maps  can  also  be  observed  in  the  direction  of  50° 
a  high  at  5.69  seconds  with  a  phase  value  of  -153°.  On 
the  contour  map  for  crosspower  between  the  horizontal  com¬ 
ponents  there  is  a  high  at  6.56  seconds  (phase  176°)  and 
in  the  coherency  map  there  is  a  high  at  6.74  seconds  (phase 
-14  7  )  ,  at  directions  of  25°  and  40°.  There  are  two  more 
nighs  on  these  maps  at  5.22  seconds  (directions  centred  at 
25  )  and  4.65  seconds  (without  a  well  defined  direction) , 
their  phases  being  -22°  and  170°  respectively. 

At  Penticton  there  is  a  peak  on  the  autopower 

spectrum  of  the  vertical  component  at  6.10  seconds.  There 

is  a  marked  high  at  this  period  on  the  contour  map  for  the 

crosspower  between  the  vertical  and  radial  components 

(between  120°  and  160°)  centred  at  140°  with  a  corresponding 

high  on  the  coherency  map  at  150°.  This  direction  actually 

falls  in  the  northeast  quadrant  at  50°  east  of  north.  Since 

there  is  a  180°  uncertainty  in  direction  using  this  method 

this  corresponds  to  a  source  in  the  Pacific  40°  south  of 

west.  This  choice  of  directions  is  also  indicated  by  the 

.  o 

low  corresponding  to  this  high  at  50  west  of  north,  which 
puts  a  low  in  the  Pacific  implying  the  absence  of  a  source 
in  that  direction  which  is  not  true.  The  phase  for  this 
event  is  50°.  A  second  high  at  6.24  seconds  due  west  can 
be  seen  on  both  crosspower  and  coherency  maps.  Phase  for 
this  event  is  36.5°.  This  high  is  actually  merging  into 
the  first  one,  and  thus  the  direction  of  approach  is  a  cone 
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the  sides  of  which  make  an  angle  of  approximately  50°.  A 
high  at  4.64  seconds  in  the  direction  of  45^  south  of  west 
can  also  be  observed  on  the  crosspower  map.  Even  though 
there  is  a  high  of  magnitude  of  0.6  on  the  coherency  map  at 
this  period,  there  is  no  well-defined  direction  for  it.  Phase 
for  this  event  is  -10°.  There  is  a  high  at  6.24  seconds  on 
the  contour  map  for  the  crosspower  between  the  two  horizontal 
components  whose  direction  is  20°  south  of  west.  There  is 
a  corresponding  high  on  the  contour  map  of  coherency  between 
these  two  components.  The  phase  for  this  is  -177°.  A  sec¬ 
ond  high  at  4.57  seconds  can  be  seen  on  both  crosspower  and 
coherency  maps  for  the  horizontal  components,  with  a  phase 
value  of  175°. 

On  the  record  chosen  for  Edmonton  there  are  two 
peaks  at  8.00  seconds  and  4.27  seconds  on  the  autopower 
spectrum  of  vertical  component.  Two  highs  (Directions  of 
which  are  25°  north  of  west  and  10°  north  of  west)  can  be 
observed  on  the  contour  map  for  the  crosspower  between  the 
vertical  and  radial  components.  Results  on  the  coherency 
map  are  very  ambiguous  and  no  clear  cut  highs  can  be  seen 
on  this  map.  Both  these  events  have  the  same  coherency  of 
0.6  with  directions  of  30°  north  of  west,  and  west.  Phases 
for  these  two  events  are  -6°  and  -11°.  A  high  at  5.33  sec¬ 
onds  can  be  seen  on  the  contour  map  for  the  crosspower  bet¬ 
ween  the  two  horizontals  at  a  direction  of  40°  north  of 
west.  Two  highs  on  the  coherency  map  for  these  two  components 
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can  be  seen  at  6.40  seconds  and  4.00  seconds  at  the  directions 
of  45  and  55°  north  of  west.  The  phases  for  these  two  are 
-168°  and  -175°. 

On  the  record  chosen  for  Resolute  thare  are  two  peaks 
on  the  autopower  spectrum  of  the  vertical  component  at  6.92 
seconds  and  5.45  seconds.  There  are  two  highs  at  6.92  and 
5.57  seconds  with  directions  of  30°  south  of  west  and  5° 
south  of  west  on  the  coherency  contour  map  for  vertical  and 
radial  components.  The  corresponding  crosspower  map  shows 
highs  at  6.92  seconds  and  5.45  seconds,  at  directions  of  15° 
north  of  west  and  30°  north  of  west.  There  is  thus  no  agree¬ 
ment  between  the  directions  obtained  from  crosspower  peaks 
and  from  coherency  peaks.  The  phases  for  these  two  events 
are  -88°  and  -102^,  indicating  Rayleigh  motion.  Three  highs 
at  6.91,  5.57  and  4.92  seconds  can  be  seen  on  the  contour 
map  of  crosspower  between  the  two  horizontals  at  directions 
of  55°  north  of  west,  40°  north  of  west  and  55°  north  of 
west.  Three  corresponding  highs  at  6.92,  5.69  and  5.02 
seconds  appear  on  the  contour  map  of  coherency  between  the 
two  horizontals  at  directions  of  55°  north  of  west,  40° 
north  of  west  and  55°  north  of  west  respectively.  Their 
phases  are  172°,  173°  and  173°  respectively.  Though 
Rayleigh  and  Love  waves  are  indicated,  their  direction  of 
approach  is  ambiguous  at  Resolute,  but  a  cone  of  direction 
of  55°  north  of  west  to  35°  south  of  west  seems  to  indicate 
a  source  in  the  Pacific  ocean.  This  large  angular  uncertainty 
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in  direction  is  also  verified  by  the  broadness  of  the 
high  in  vertical-radial  crosspower. 

The  autopower  of  the  vertical  component  peaks  at 
6.40  seconds  at  Ottawa.  The  high  on  the  contour  map  of 
the  crosspower  between  the  vertical  and  radial  components 
occurs  at  60°  north  of  west.  A  high  on  the  coherency  map 
for  these  two  components  appears  at  6.56  seconds  with  a 
direction  of  45°  north  of  west.  The  phase  for  this  event 
is  -179°.  There  are  four  highs  on  the  contour  map  of  cross¬ 
power  between  two  horizontal  components  at  periods  of  7.11, 
6.24,  5.57  and  4.92  seconds.  Their  directions  are  80°  north 
of  west,  70°  north  of  west,  55°  north  of  west  and  65°  north 
of  west.  Three  of  these  are  supported  by  highs  on  the  coh¬ 
erency  map  for  the  two  horizontals  occurring  at  7.31,  5.69 
and  4.92  seconds.  Their  directions  are  north,  60°  north  of 
west  and  also  60°  north  of  west,  respectively.  In  addition 
there  is  a  plateau  of  0.3  coherency  between  6.74  -  6.24 
seconds  with  no  well-defined  direction.  The  phases  for  all 
these  events  are  close  to  -180°  or  0  . 

Vi .  4  The  contour  diagrams  in  general  show  highs  in  the 

North  West  or  South  West  quadrants  indicating  a  source  in 
the  Pacific  Ocean.  Lows  are  well  marked  (except  in  one  or 
two  cases)  and  90°  away  from  the  highs  in  the  case  of  vertical- 
eadial  system  and  45°  in  the  case  of  radial-transverse  sys¬ 
tem  as  expected.  In  certain  cases  lows  are  more  clearly 
defined  than  the  highs,  which  tend  to  be  broad. 


; 
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A  comparison  of  the  four  maps  at  a  given  location 
shows  a  general  agreement  between  the  directions  of  approach. 
It  is  only  at  Resolute  that  the  directions  obtained  from 
the  crosspower  highs  and  the  coherency  highs  for  the  vertical- 
radial  system  differ  by  as  much  as  40°,  otherwise  they  are 
within  +10  range.  Not  only  do  the  directions  obtained  from 
two  quantities  in  one  set  agree  to  +  10°  but  they  are  also 
close  to  those  obtained  from  the  other  set,  viz.  the  radial- 
transverse  system.  It  can  be  said  therefore  that  the  direc¬ 
tions  of  approach  of  Rayleigh  and  Love  waves  differ  at  most 
by  10°.  The  discrepancy  at  Resolute  is  not  readily  explain¬ 
able.  As  will  be  shown  below  the  directions  obtained  here 
can  be  favourably  compared  with  those  of  the  storm  locations. 

Besides  the  directions,  magnitudes  of  highs  and  lows 
offer  very  interesting  features.  It  was  pointed  out  earlier 
that  the  ratio  between  high  and  low  provides  a  measure  of 
the  signal  to  noise  ratio  in  the  microseisms.  Two  such  ratios 
exist  in  the  two  sets  of  computations;  one  corresponds  to 
the  crosspower  and  the  other  to  coherency.  The  crosspower 
ratios  are  as  low  as  1.33  (for  the  7.31  second  event  at 
Ottawa)  and  as  high  as  22  (for  the  5.45  second  event  at 
Resolute) ;  the  average  figure  however  seems  to  be  around  8. 
Thus  the  amplitude  signal  to  noise  ratio  is  about  3.  The 
range  between  the  signal  and  noise  level  therefore  varies 
from  as  low  as  3db  to  well  over  12db.  This  comes  as  no 
surprise  in  view  of  the  sinusoidal  appearance  of  the  micro- 
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seisms  on  the  records. 

In  most  cases  the  noise  level  for  coherency 
is  about  0.1/  although  in  one  case  it  is  as  large  as  0.5. 

But  the  very  fact  that  the  coherency  in  the  direction  of  the 
arriving  waves  rarely  exceeds  0.6  for  the  vertical-radial 
system  coupled  with  the  occurence  of  spurious  phase  values, 
shows  that  other  types  of  waves  are  recorded,  especially 
on  the  vertical  component.  (Note  that  the  coherency  values 
for  radial- transverse  system  are  usually  greater  than  0.8). 
This  conclusion  is  reinforced  by  the  vertical  to  radial 
ratios  calculated  from  these  diagrams  which  rarely  agree 
with  theoretical  ratios. 

There  are  three  cases  where  a  Rayleigh  type  of 
motion  could  be  identified.  Two  of  them  are  at  Resolute 
where  the  phase  values  indicate  a  retrograde  motion;  the 
third  at  Victoria  where  the  phase  value  indicates  a  prograde 
motion.  At  Resolute,  for  this  sample  record,  there  is  an 
excellent  agreement  between  the  observed  and  theoretical 
values  of  the  vertical  to  radial  ratios  supporting  the 
identification  of  a  Rayleigh  wave.  At  the  other  stations 
(including  Victoria)  the  vertical  component  is  usually 
larger  than  expected,  and  the  phase  is  seldom  appropriate 
for  Rayleigh  motion.  If  one  were  to  interpret  all  these 
features  (large  values  for  vertical  amplitude,  wrong  phase 
values,  and  yet  high  signal  to  noise  ratio)  as  due  to  the 
presence  of  long  period  body  waves  (especially  leaking  P- 
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TABLE  VI. la 


Period  Components 


Direction 


(Secs)  V-Vertical  Power 
R-Radial 

T-Transverse  ~  " 
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TABLE  VI. lb 


Period  Components 


Magnitude  Station 


(Secs) 

V-Vertical 

Power 

Coherency 
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Coherency 

R- Radial 

T- 

•Transverse 
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Low 
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6.10 

V 

60 

hd 

M 

2 

6.24 
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H 
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TABLE  VI. lc 


Period 

Components 

V/R 

V/R  Station 

(Secs) 

V-Vertical 

R- Radial 
T-Transverse 

Observed 

Theoretical 

6.24 

• 

V-R 

2.9 
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0.72  OTTAWA 
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modes) ,  one  has  to  accept  that  the  record  at  Resolute  is 
unusual,  and  thus  that  pure  Rayleigh  waves  in  microseisms 
are  the  exception  rather  than  the  rule. 

In  table  VI. -1  are  summarized  the  several  quantities 
computed  from  tnese  diagrams  and  discussed  above  in  sections 
VI. 3  and  VI. 4. 


VI .  5  In  summary  then  the  following  features  can  be 
recognized  from  these  contour  diagrams: 

1.  As  a  rule  the  contours  tend  to  be  broad 
in  all  cases.  There  are  cases  when  the 
contours  do  centre  around  a  small  closed 
area  but  these  are  only  exceptions. 

2.  Contour  lines  are  broader  for  the  radial- 
vertical  system  than  for  the  radial- 
transverse  system. 

3.  The  broadness  of  the  contours  is  indicative 
of  uncertainty  in  direction  of  approach. 

The  farther  the  station  is  removed  from 
the  coast,  the  broader  are  the  contours. 

4.  There  is  no  marked  improvement  in  the 
values  of  coherency  even  after  rotation. 

5.  Cross  sections  of  the  contours  on  the 
cross  power  grids  do  show  very  well  marked 
peaks,  a  deduction  made  from  the  closely 
bunched  nature  of  the  contours. 

6.  The  lows,  the  highs  and  the  plateaux  are 
very  well  shown  on  these  diagrams. 

7.  Highs  and  lows  are  separated  by  90°  on 
radial-vertical  contour  maps  and  are  re¬ 
peated  every  full  rotation  (180°) .  The 
highs  and  lows  on  the  contour  maps  for 
two  horizontal  components  are  separated 
by  45°  and  they  are  repeated  every  half 
rotation  (90°) . 


* 
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8.  A  comparison  was  made  between  the  direc¬ 
tions  obtained  from  rotated  coherency 
and  rotated  cross  power  with  the  direc¬ 
tions  obtained  from  the  bearings  of  storm 
locations.  These  comparisons  are  shown 
in  Table  VI . 2 .  The  uncertainty  in  the 
direction  increases  as  one  moves  away 
from  the  coast.  The  best  fit  is  for 
Victoria,  the  next  best  for  Penticton,  in 
that  order  Edmonton,  Resolute  and  the 
worst  for  Ottawa.  The  discrepancy  ranges 
from  as  low  as  -3  to  as  high  as  54°, 
but  considering  the  limitations  on  the 
data,  one  has  to  admit  that  there  is  a 
good  agreement  between  the  directions  ob¬ 
tained  from  cross  power  spectra  and 
coherency. 

9.  From  3  and  8  above,  it  can  be  deduced 
that  the  seismic  waves  excited  in  the 
ocean  bottom  by  the  overlying  water  (due 
to  second  order  pressure  fluctuations) 
are  channeled  in  the  oceanic  crust,  which 
acts  as  a  wave  guide,  until  the  energy 
reaches  the  coast  at  different  points  and 
is  then  propagated  in  the  continental 
crust  along  several  different  paths. 

This  effect  of  the  coastline  becomes 
more  pronounced  at  sites  very  remote  from 
the  coast  indicating  that  the  coastline 
acts  as  a  major  source  of  waves.  Thus 
the  simple  Longuet-Higgins  theory  needs 

to  be  modified  to  take  into  account  inter¬ 
actions  of  the  primary  Rayleigh  waves 
with  the  coastline  and  the  subsequent 
generation  of  other  wave  types. 

10.  From  the  phase  values,  from  the  signal 
and  noise  levels  as  depicted  by  coherency 
values  and  from  the  power  recorded  in 
the  vertical  component,  it  is  evident 
that  body  waves  (P-waves  in  particular) 
of  around  6-7  seconds  are  present  in 

the  microseisms. 

11.  Evidences  for  the  presence  of  surface 
waves  are  cited  where  ever  possible 
but  they  do  not  seem  to  form  a  dominant 
component  of  the  microseisms. 
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12.  The  estimation  of  proportion  of  body 
waves  in  microseisms  as  compared  to 
the  surface  waves  is  not  feasible  with 
the  information  available  at  present. 

An  array  of  seismic  stations  separated 
by  distances  comparable  to  the  wave-* 
length  of  the  surface  waves  would  yield 
phase  velocities  of  different  waves  which 
could  then  be  used  as  criteria  for 
this  estimation. 


CHAPTER  VII 


GENERAL  DISCUSSION 

VII.  1  An  attempt  thus  has  been  made  to  study  microseisms 
recorded  at  five  seismological  observatories  in  Canada  on 
dates  when  an  atmospheric  low  pressure  area  was  present 
on  the  Pacific  Ocean.  Spectral  estimates  were  mainly 
utilized  in  identifying  the  wave  types  present  in  micro¬ 
seisms,  in  determining  their  direction  of  approach  and  in 
understanding  the  mechanism  of  the  propagation  over  long 
distances.  The  data  were  Butterworth  band  pass  filtered  to 
draw  Jensen's  diagrams  and  particle  motion  trajectories. 

No  assumptions  were  made  about  the  nature  of  the 
microseisms  and  the  problem  was  approached  with  an  entirely 
open  mind.  The  results  are  revealing  and  at  times  dis¬ 
couraging:  discouraging  because  of  the  limitations  on  the 

data,  and  not  because  of  any  method  employed  in  under¬ 
standing  them.  The  data  do  not  seem  to  be  suited  for  con¬ 
structing  the  Jensen's  diagrams  because  of  the  very  com¬ 
plex  nature  of  the  microseisms  studied  here.  That  this 
method  has  failed  can  be  construed  as  an  evidence  of  micro¬ 
seisms  not  containing  solely  surface  waves.  The  presence 
of  body  waves  (especially  P-modes)  in  the  period  range  of 
6  —  7  secs  and  of  4-5  secs  was  suggested  on  the  basis  of 
several  computations.  The  evidences  for  the  presence  of 
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body  waves  in  inicroseisins  can  be  cited  in  their  order  of 
occurrence  in  this  work  as  follows: 

1.  Sequence  of  peaks  in  the  velocitv  sensitivity 
spectra  indicating  a  resonance  phenomenon 
associated  with  the  body  waves  (Phinney,  1964) . 

2.  Phase  values  between  the  vertical  and  radial 
components  which  are  usually  0°  or  180°. 

3.  Amplitude  ratios  of  vertical  to  radial  com¬ 
ponents  which  are  very  high  compared  to  the 
theoretical  values  suggesting  that  in  addition  to 
Rayleigh  waves,  long  period  P  modes  are  recorded 
on  the  vertical  component. 

4.  Jensen's  diagrams  which  show  a  tremendous 
amount  of  scattering  for  the  velocity  vectors 
pointing  that  no  simple  Rayleigh  and  Love 
waves  are  being  recorded  on  these  instruments. 

On  the  other  hand  several  unmistakable  surface 
wave  examples  were  mentioned.  It  was  also  suggested  that 
in  addition  to  fundamental  modes,  higher  modes  of  both 
Rayleigh  and  Love  waves  were  present  in  microseisms  super¬ 
imposed  in  a  complex  fashion.  In  general  surface  waves  have 
two  dominant  periods  one  between  6. 0-7.0  sec  and  the  other 
around  4.5  sec.  A  few  examples  of  2.5  second  surface 
waves  could  be  identified. 

VII. 2  Cross-power  and  coherency,  as  functions  of  two 
variables,  viz.  period  and  azimuth,  were  computed  and 
plotted  on  a  grid,  using  the  resolved  horizontal  traces. 

This  grid  was  then  used  to  construct  contour  diagrams.  Ihe 
purpose  of  this  was  to  determine  the  direction  of  approach 
of  microseismic  waves.  The  results  were  gratifying  in 
so  far  as  it  was  found  that  direction  Oi  approach 
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could  be  determined  to  an  accuracy  of±10°.  Comparison 
of  the  directions  obtained  from  cross  power  and  coherency 
between  the  vertical  and  horizontal  components  provides  a 
quadruple  check  on  this  measurement.  Not  always  did  these 
four  directions  agree,  but  seldom  did  they  have  a  dis¬ 
crepancy  of  more  than  10°.  These  directions  were  then  com¬ 
pared  with  those  obtained  from  the  storm  locations.  Con¬ 
sidering  the  limitations  on  the  data,  the  agreement  seems 
to  be  good. 

VII. 3  It  is  suggested  here  that  the  seismic  waves  excited 
in  the  ocean  bottom  by  the  second  order  pressure  fluctua¬ 
tions  in  the  overlying  water  (Longuit-Higgins ,  op.  cit.) 
are  propagated  in  multiple  directions  in  the  oceanic  crust, 
which  acts  as  a  wave  guide.  The  energy  is  then  transferred 
from  the  oceanic  crust  to  the  continental  crust  and  in  so 
doing  undergoes  an  interaction  with  the  coastline  which 
possibly  acts  as  a  number  of  secondary  sources  from  which 
the  waves  travel  in  the  continent,  so  that  a  recording 
site  would  look  at  different  waves  approaching  the  site 
in  different  directions.  The  farther  is  the  site  from  the 
coast,  the  wider  are  the  angles  and  hence  the  more  uncertain 
is  the  direction  of  approach.  It  was  found  that  the  ver¬ 
tical  component  of  the  microseismic  waves  attenuates  very 
nearly  as  —  and  the  horizontal  approximately  as  —  also, 
where  r  is  the  distance  from  the  source. 
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